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INTRODUCTION 

The  widespread  use  of  square  hysteresis  loop  ferrite  materials  in 
computer  technology  has  directed  attention  to  the  problem  of  constructing  an 
engineering  model  of  the  process  of  flux  reversal.   Such  a  model  would, 
ideally,  allow  the  square  loop  ferrite  specimen  to  be  treated,  in  circuit 
design,  as  any  other  common  circuit  component.   For  the  simplest  geometry, 
the  toroid,  the  problem  is  just  to  establish  the  voltage-current  relation 
for  a  two-terminal  network,  the  terminals  being  the  two  ends  of  a  conductor 
which  links  the  magnetic  flux  of  the  toroid. 

Notable  success  has  been  achieved  by  previous  investigators.   A 
broad  area  of  theoretical  and  experimental  work  was  effectively  coordinated 
by  Menyuk  and  Goodenough  (1955)  resulting  in  basic  equations  of  the  process 
of  flux  reversal  by  domain-wall  motion  and  clarifying  certain  questions  of 
the  related  microscopic  details.   First,  Haynes  (1958)  and,  later,  Lindsey 
(1959)  specified  the  function,  F(<r>),  which  remained  in  the  Menyuk- Goodenough 
equations,  resulting  in  practical  equations  of  flux  reversal. 

The  simplest  experiment  producing  flux  reversal  in  a  magnetic  core 
consists  of  the  application  of  a  step  function  of  mmf  with  the  core 
initially  in  one  of  its  two  reproducible  (saturated)  states.   Even  for  this 
experiment,  the  known  models  are  demonstrably  inaccurate,  although  adequate 
for  some  applications.   It  will  be  shown  presently  that  this  problem  may 
be  resolved  by  generalizing  the  models  of  Lindsey  and  Haynes. 

The  more  interesting  problem  which  is  investigated  here  is  that 
of  construction  of  a  model  for  "partial  switching",  that  is,  the  case  of 
flux  reversal  proceeding  from  an  initially  unsaturated  state.   Attention  is 
focused  on  a  particular  type  of  partial  switching,  that  in  which  the  inter- 
mediate state  is  reached  by  a  pulse  of  mmf  in  the  "UP"  (set)  direction. 
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The  response  to  a  pulse  in  the  "DOWN"  (reset)  direction  is  then  examined. 
Certain  invariant  and  semi -invariant  quantities  are  found  to  exist  for  the 
observed  response.   In  particular,  the  envelope  of  peak  amplitudes  of  the 
voltage  versus  time  curves  itself  reaches  a  maximum,  the  "maxpeak",  which 
appears  at  first  sight  to  he  anomalously  large.   It  is  demonstrated  that  a 
logical  extension  of  Lindsey's  model  is  adequate  to  explain  the  observations 
in  a  rudimentary  way.   In  addition,  the  study  of  the  maxpeak  effect  model 
is  believed  to  provide  a  more  stringent  test  of  some  of  the  assumptions 
of  the  earlier  models . 

In  the  construction  of  flux  reversal  models,  the  principal  con- 
straint in  the  present  work  is  that  the  model  be  based  on  microscopic 
switching  processes  which  are  physically  rational.   Stated  negatively, 
purely  mathematical  models  are  not  considered. 


I.   FUNDAMENTAL  PROPERTIES  OF  MAGNETIC  MATERIALS 
1.1  Energy  Relations 

In  order  to  proceed  with  the  discussion  of  magnetic  effects,  it  is 
desirable  to  review,  however  briefly,  the  kinds  of  energy  upon  which  these 
effects  depend.   This  has  been  done  by  Van  Vleck  (1956),  whose  classification 
is  followed  here . 

(1)  Zeeman  Energy 

This  is  the  energy  which  exists  in  the  interaction  of  the 
applied  magnetic  field,  H  ,  with  electron  spin  and  orbital 
angular  momenta  which  can  be  written 

Ez   =   -gp  E  (H.  S.)  (l.i) 

i 

where  g  is  the  spectroscopic  splitting  factor  (  ~  2) 

p  is  the  Bohr  magneton 

■*"  th 

S.  is  the  spin  angular  momentum  of  the  i   electron. 

Since  the  electrons  tend  to  occur  in  antiparallel  pairs  in 

conformance  with  the  exclusion  principle,  only  a  very  small 

percentage  of  electrons  in  a  material  absorb  net  energy  from 

the  applied  field.   The  substantial  magnetic  effects  in  Fe, 

Ni,  and  Co  are  due  to  the  unusually  large  number  (up  to  5) 

of  unpaired  electron  spins. 

(2)  Classical  Electromagnetic  Dipolar  Energy 

This  is  the  energy  effect  commonly  observed,  on  macro- 
scopic level,  as  the  tendency  of  bar  magnets  to  align 
themselves  such  that  the  direction  of  maximum  length  corres- 
ponds to  the  direction  of  the  magnetization  vector.   At 
atomic  distances,  the  dipolar  energy  is  masked  by  other 
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interactions.   Its  long  range  effect,  however,  accounts  for 
the  appearance  of  domain  structure. 
(3)  Exchange  Energy 

The  exchange  energy  is  most  simply  expressed  by  the 
relation 

E   =   -2J   Z    S.  •   S.  (1.2) 

where  J  is  the  quantum  mechanical  exchange  integral  and 
i'**j  implies  summation  over  nearest  neighboring  atoms. 
This  energy  accounts  for  a  large  portion  of  the  short-range 
ordering  of  spins  which,  it  should  be  noted,  is  a  prime 
requisite  to  establishing  magnetic  effects.   It  can  also  be 
seen  that  exchange  energy  is  independent  of  external 
direction  references,  i.e.,  it  depends  solely  upon  the 
relative  spin  directions  and  magnitudes. 
(h)     Anisotropic  Exchange  Energy 

This  energy  results  from  the  interaction  of  spin  with 
spherically  asymmetric  orbital  charge  distribution.   This 
effect  is  similar  to  the  classical  dipolar  energy  except 
that  anisotropic  exchange  is  sufficiently  stronger  at 
small  distances  to  account  for  the  powerful  anisotropic 
effects  observed  in  solids.   The  anisotropy  energy  has 
the  macroscopic  form: 

EA  =  Ax  +  Bc^  +  Cc£  (1.3) 

where  A,  B,  and  C  are  constants  and  a   ,  a   ,  and  a  are  the 
direction  cosines  of  the  angles  between  the  direction  of 
magnetization  and  the  principal  crystal  axes. 


A  sketchy  physical  conception  of  the  origin  of  magnetic 
effects  may  be  extracted  from  the  more  formal  statements 
above.   In  a  magnetic  specimen,  the  crystal  structure  deter- 
mines the  favorable  orbital  electron  motion.  The  orbital 
motion  affects,  by  anisotropic  exchange,  the  spin  direction. 
In  ferromagnetic  materials,  a  substantial  net  spin  may 
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Figure  1.1 


A  sketchy  physical  conception  of  the  origin  of  magnetic 
effects  may  be  extracted  from  the  more  formal  statements 
above.   In  a  magnetic  specimen,  the  crystal  structure  deter- 
mines the  favorable  orbital  electron  motion.  The  orbital 
motion  affects,  by  anisotropic  exchange,  the  spin  direction. 
In  ferromagnetic  materials,  a  substantial  net  spin  may 
result . 

1.2  Domains  and  Bloch  Walls 

In  a  perfect,  finite,  single  crystal  of  a  magnetic  material,  it 
might  appear  that  a  net  magnetization  exists  which  would  necessarily 
generate  return  flux  paths  through  the  surrounding  medium,  e.g.,  air. 
Due  to  the  relatively  large  energy  required  to  establish  the  flux  in  a 
low  permeability  medium,  it  is  preferable  to  close  the  field  within  the 
magnetic  medium.  This  can  be  done,  for  example,  as  illustrated  in  Fig.  1.1 
a,  b  where  the  magnetization  vectors  are  indicated.   In  other  terms, 
the  surface  pole  distribution  (demagnetization  component)  due  to  the 
boundary,  affects  the  interior  by  the  long-range  dipolar  mechanism. 
The  resulting  regions  of  uniform  magnetization  are  domains  or,  "Weiss 
domains".  The  boundary  between  domains  is  a  Bloch  wall  (domain  wall) 
whose  characteristics  are  explained  in  the  following  paragraph. 
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Figure  1.1 
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If  the  maximum  space  rate  of  change  of  the  direction  of  magnetiza- 
tion in  the  Bloch  wall  is  very  great,  it  can  he  seen  from  equation  (1.2)  that 
the  exchange  energy  would  he  large  (J  >  0  assumed)  so  that  this  energy  favors 
a  wide  Bloch  wall  with  negligible  magnetization  change  from  one  atom  to  the 
next.  However,  assuming  only  the  +x  and  +y  directions  of  Fig.  1.1  correspond 
to  low  anisotropy  energy,  the  slow  change  of  magnetization  in  the  Bloch  wall 
would  compel  large  numbers  of  atomic  magnetization  vectors  into  directions 
of  high  anistropy  energy.  Thus  a  compromise  arises  which  results  in  Bloch 
walls  having  small  hut  non-negligible  widths  (of  the  order  of  100  A). 

The  change  in  angle  of  the  magnetization  vector  from  one  domain 
to  another  contiguous  domain  is  used  to  classify  domain  walls.   In  Fig.  1.1  a, 
for  example,  there  are  four  90  and  one  180  domain  walls . 

The  characteristics  of  domains  will  be  relegated  to  section  II. 1. 

1.3  Ferrites 

The  ferrites  are  metallic  oxides  whose  molecular  formula  is 
X0  •  Y  0  where  X  corresponds  to  a  metallic  atom  of  valence  2  and  Y  corresponds 
to  a  metallic  atom  of  valence  3. 

A  wide  variety  of  metals  may  be  employed,  the  iron  group  being 
particularly  common.  The  oxygen  atoms  provide  bonds  for  the  metallic 
electrons  which  would  otherwise  be  free  for  electrical  conduction.  Thus, 
the  electrical  conductivity  of  ferrites  may  be  very  low.  This  fact  is  of 
engineering  importance  in  two  ways.   First,  the  eddy  current  losses  which 
characterize  other  magnetic  materials  are  eliminated.   Secondly,  the  maxi- 
mum rate  of  change  of  magnetization,  which  would  otherwise  be  limited  by 
the  eddy  current  field,  may  be  very  great.  The  second  factor  is  particularly 
desirable  in  computer  technology. 

In  ferromagnetic  materials,  all  unpaired  electron  spins  within 
a  unit  cell  of  the  material  are  approximately  aligned  so  that  all  contribute 
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to  the  magnetization.   In  ferromagnetic  materials,  (i.e.,  ferrites),  certain 
electron  spins  which  are  unpaired  within  their  own  atom,  may  cancel  the 
magnetic  contribution  of  similar  electrons  in  neighboring  atoms.   If  the 
cancellation  within  a  unit  cell  is  total,  the  material  is  said  to  be  anti- 
ferromagnetic.   If  the  cancellation  is  only  partial,  the  material  displays 
ordinary  ferrimagnetism. 

The  metallic  atoms  which  can  contribute  to  the  magnetization  in 
ferrites  may  occupy  two  types  of  sites,  called  A  and  B  sites.   In  a  unit 
cell  of  ferrite,  the  number  of  such  sites  occur  in  the  ratio  1:2.   If  the 
exchange  energy  between  B  sites  is  large,  a  degree  of  antiparallel  align- 
ment of  B  atoms  results.   Thus,  the  attendant  magnetic  contribution  of 
the  B  sites  is  small.   This  structure  is  called  "normal".   If  the  B  atomic 
spins  are  aligned  but  antiparallel  to  the  spin  of  A  atoms,  the  structure 
is  said  to  be  of  the  "inverse"  type. 

1.4  Preparation  of  Ferrites 

The  ferrite  bodies  which  are  of  interest  here  have  polycrystal- 
line  structure  and  are  prepared  by  the  following  procedure.   A  mixture  of 
metallic  oxides  is  powdered  by  a  ball-milling  process;  the  resulting 
mixture,  including  possibly  a  binder  material,  is  packed  into  a  mould  and 
the  desired  hole(s)  punched  out.   The  body  is  then  sintered  for  periods  of 
one  minute  to  several  hours  at  temperatures  of  approximately  1300  C  in  the 
presence  of  a  gaseous  atmosphere,  such  as  nitrogen.  The  various  parameters 
must,  in  general,  be  rigidly  controlled  in  order  to  obtain  adequately 
reproducible  results. 

It  should  be  noted  that  no  magnetic  orientation  (annealing)  is 
employed  so  that  the  bodies  are  magnetically  isotropic  barring  possible 
shape  effects. 
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As  a.   result  of  its  method  of  fabrication,  a  ferrite  body  usually 
contains  a  number  of  voids  (typically  10$  of  the  volume)  whose  distribution 
in  shape,  size,  and  position  appears  to  be  quite  random.  The  crystal 
structure  is  further  disarrayed  by  grain  boundaries  which  leave  a  random 
distribution  of  volumes  having  a  periodic  crystal  structure. 

Photomicrograms  of  the  cross-section  of  ferrite  bodies  are  found 
in  Smit  and  Wijn  (1959). 


II.   FLUX  REVERSAL 

II. 1  Engineering  Model  of  Flux  Reversal 

In  order  to  precisely  design  electrical  circuits  using  a  toroidal 
magnetic  core,  it  is  only  necessary,  theoretically,  to  know  the  current- 
voltage  relation  for  a  single  wire  linking  the  toroidal  path  for  all 
possible  (time -varying)  currents  through  the  wire.   From  the  viewpoint  of 
practical  engineering  design,  such  a  relation  would  probably  be  too 
cumbersome  for  general  use.   It  would,  however,  be  desirable  to  have  a 
mathematical  model  of  flux  switching  in  these  cores  which  is  reasonably 
accurate  over  a  wide  range  of  conditions,  even  at  the  expense  of  being 
too  cumbersome  for  usual  design  work.   The  mathematical  model  would  then 
have  the  possibility  of  being  simplified  for  special  cases. 

Previous  efforts  to  formulate  such  models  have  been  based  on 
the  assumption  that  the  magnetic  specimen  is  initially  in  a  saturated 
state.   Such  equations  are  of  the  form: 

M  =  f(M,H) 

with     M  =  the  time  rate  of  change  of  magnetization 

M  =  the  magnetization,  and 

H  =  the  applied  field  intensity. 
This  equation  can  also  be  written 

0   =   f(0,  H) 


to  a  negligible  approximation  in  most  cases  since  the  induction  directly 
contributed  by  H,  (\x    H  ),  is  very  small  compared  to  M  in  the  relation 

B  =  u  H  +  M. 
o 
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Only  at  times  when  H  is  large  is  there  any  practical  concern  in  replacing  M 
by  0. 

The  circularly  directed  component  of  magnetization  is  the  only 
component  of  interest  here;  hence,  scalar  notation  is  employed. 

The  equations  of  the  previous  type  which  have  been  derived  are 
based  on  switching  by  (a)  180  domain-wall  motion,  (b)  rotation.   The  equation 
for  rotation  was  derived  by  Gyorgy  (1957)  from  the  Landau-Lifshitz  phenomeno- 
logical  damping  equation.   The  present  investigation  will  be  limited  to  the 
domain-wall  motion  model  formulated  by  Menyuk  and  Goodenough  (1955)  and  ex- 
tended by  Haynes  (1958)  and  Lindsey  (1959). 

II. 2  Rotation 

Flux  changes  may  occur  by  rotation  as  a  result  of  the  interaction  of 
the  magnetic  field  intensity,  H,  with  the  magnetization,  M.  The  resulting 
torque,  causing  rotation,  is  T  =  M  x  H. 

If  the  field,  H,  is  small,  only  individual  atomic  spins  are  rotated 
against  the  anisotropic  exchange  force  and  the  action  is  generally  characterized 
by  small,  fast  ( &   1  mas)  reversible  changes  of  magnetization.   If,  however, 
sufficiently  large  fields  are  applied,  the  magnetization  of  a  complete  domain 
or  sample  of  material  may  be  moved  past  the  direction  of  maximum  anisotropic 
or  demagnetizing  force  such  that,  upon  removal  of  the  applied  field,  a 
different  arrangement  of  domains  exists,  i.e.,  the  flux  change  is  now 
irreversible.   At  some  intermediate  field,  it  is  possible,  according  to 
Gyorgy  (1958),  to  rotate  the  magnetization  in  such  a  way  that  the  demagnetizing 
effects  may  be  neglected.   If  the  magnetization  of  the  whole  volume  of  material 
in  question  is  irreversibly  changed  by  rotation  against  demagnetizing  forces, 
the  switching  process  is  said  to  be  coherent  or  uniform  rotational  switching. 
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If  the  magnetization  is  rotated  in  a  way  such  that  the  demagnetizing  effects 
may  be  neglected,  the  switching  process  is  then  incoherent  or  nonuniform 
rotation. 

Estimates  by  Smith  (i960)  indicate  that  the  minimum  drive  for  which 
rotational  effects  can  occur  is  four  times  the  threshold  field  for  materials 
similar  to  the  ones  used  in  the  present  study.   The  present  investigation 
will  be  carried  out  for  fields  smaller  than  the  rotation  threshold. 

II. 3  Domain  Wall  Motion 

The  voids,  grain  boundaries,  and  other  imperfections  in  a  poly- 
crystalline  ferrite  have  a  fundamental  effect  on  the  magnetic  properties  of 
the  material.   It  can  readily  be  seen  that  voids  of  sufficient  size  are 
analogous  to  the  boundary  of  the  material,  thereby  causing  the  formation  of 
90  domain  walls  (closure  domains)  and  surface  poles  as  explained  in  section 
1.2. 

The  process  of  flux  reversal  by  domain  wall  motion  is  generally 
pictured  as  follows:   starting  from  an  initial  condition  of  spontaneous 
magnetization,  a  sufficiently  large  magnetic  field  applied  in  the  reverse 
direction  causes  a  180  Bloch  wall  to  be  initiated  or  "nucleated"  along  a 
grain  boundary.   As  the  Bloch  wall  moves,  the  volume  in  its  wake  is  left 
with  the  magnetization  reversed.   The  wall  continues  to  move  until  either 

(1)  it  meets  another  wall  moving  toward  it,  which  annihilates  both,  or 

(2)  the  driving  field  is  removed,  in  which  case  the  wall  comes  to  rest  in 
a  way  consistent  with  requirements  of  minimum  potential  energy. 

As  a  domain  of  reverse  magnetization  grows  due  to  the  pressure  of 
the  applied  magnetic  field,  the  increasing  domain  wall  surface  area  absorbs 
energy  from  the  driving  field.   In  high  resistivity  material,  this  energy  is 
lost  to  damping  effects  which  arise  from  the  delayed  response  of  atomic  spins 
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to  the  force  which  tends  to  change  their  direction.   It  is  usually  assumed 
that  this  damping  is  viscous  so  that  domain  walls  subjected  to  a  steady  field 
will  move  with  constant  velocity,  once  the  initial  acceleration  has  been 
completed . 

The  conditions  which  must  be  met  in  order  to  attain  squareness  of 
the  hysteresis  loop  are  also  of  interest.   According  to  Goodenough  (195^)^ 
the  nucleation  of  a  domain  wall  occurs  when  the  energy  required  to  form  surface 
poles  at  the  grain  boundary  becomes  larger  than  the  energy  required  to  form  a 
region  of  reverse  magnetization.   Therefore,  at  applied  field  intensities 
below  the  nucleation  threshold,  the  formation  of  surface  poles  must  be  pre- 
ferred.  This  is  expressed  by  Goodenough' s  result, 

LIC (cos  9.,  -  cos  9-)  <  60  6~ 

s v     1        2         w 

where   L  is  the  mean  grain  diameter 

M  is  the  spontaneous  magnetization 

9  ,  9  are  the  directions  of  magnetization  in  adjacent  grains 

relative  to  their  common  boundary  and 

<S     is  the  surface  energy  density  of  a  180  Bloch  wall. 
w 

The  above  condition  is  believed  to  be  met  in  magnetically  annealed  materials 

because  of  the  small  value  of  (cos  9  -  cos  9  )  while  in  polycrystalline 

ferrites,  it  is  due  to  the  small  value  of  W . 

s 

II. h     Review  of  Menyuk  and  Goodenough' s  Contribution 

A  major  step  toward  the  development  of  a  model  of  flux  switching 
was  taken  by  Menyuk  and  Goodenough  (1955)-  They  consider  the  expanding 
domains  to  be  ellipsoids  having  large  semimajor  to  semiminor  axis  ratio  and 
derive  the  relation 

P  d  T~  =  21  (H  -  H  )  <cos  9>  (2.1) 

dt        s  m    o 
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where  f3  is  the  damping  parameter 

<r>  is  the  mean  value  of  semiminor  axis  at  the  center  of  the  domain 

I  is  the  spontaneous  magnetization 

H  is  the  applied  field  intensity 
m 

H  is  the  threshold  field,  that  is,  the  upper  bound  of  field 

intensity  for  which  the  domain  wall  velocity  is  zero,  and 

9  is  the  angle  between  the  domain's  magnetization  vector  and 

the  applied  magnetic  field  direction. 

It  should  be  noted  that  the  preceding  equation  is  based  on  the 

assumption  that  the  domain  walls  move  with  constant  velocity,  that  is, 

H  -  H  =  constant  > 0.  Integration  of  (2.1)  yields  the  now-classic  relation 
mo 

for  H  independent  of  time: 

(H  -  H  )Z         =     S  (2.2) 

x  m    o   s      w  v 

where  t     is  the  switching  time  and 
cs  to 

S  is  the  switching  coefficient,  a  constant  which  depends  on  the 
properties  of  the  core  material. 

They  proceed  further  to  write  the  time  rate  of  change  of  irreversible 
flux,  formally,  as 


d0 

9A 


°C  F(<r>)(H  -  H  )  (2.3) 

dt  v     x  m    o/  v  ~JJ 

irr 


where  F(<r>)  -  \-rr~     >  A  being  the  cross-sectional  area  of  the  core  in  which 
flux  has  been  reversed. 

A  typical  example  of  a  full  switching  curve  is  given  in  Fig.  2.1. 

A  direct  quotation  from  the  Menyuk-Goodenough  paper  will  serve  to- 

elucidate  their  view  of  the  flux  reversal  process . 

"F(<r>)  measures  the  rate  of  change  of  flux  as  a  function 
of  the  constant-velocity-wall  positions.   Since  the  average  wall 
velocity  is  assumed  to  accelerate  to  a  constant  velocity  in  a 
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time  t  «  t ,  F(<r>)  is  a  maximum  when  the  wall  area  encompassing 
the  growing  domains  is  a  maximum.   The  wall  area  surrounding  a 
growing  domain  increases  with  r  until  it  collides  with  similarly 
growing  domains  and  is  annihilated.   Because  of  the  ramdom 
distribution  of  the  growing  domains  throughout  the  specimen,  at 
larger  values  of  r,  F(<r>)  decreases  to  zero  with  a  Gaussian-like 
tail." 


II. 5   Haynes '  Model 

Haynes  (1958),  proceeding  from  the  previously  described  development 
of  Menyuk  and  Goodenough,  introduced  assumptions  required  to  specify  the  function, 
F(<r>).   Haynes  states  the  following  explicit  assumptions: 

(1)  Only  irreversible  wall  motion  is  treated. 

(2)  The  nucleation  field  strength  is  equal  to  the  threshold  field. 

(3)  The  number  of  nucleated  domains  is  independent  of  H  -  H  . 
(h)     The  nucleation  process  and  the  finite  minimum  size  of  a 

domain  are  neglected. 

(5)  The  variation  of  H  with  radial  position  with  a  core  is  neglected. 

(6)  3A  /3<s>  is  a  function  of  domain  wall  position,  <s>,  and  not  of 
the  velocity  d<s>/dt. 

Haynes  makes  the  further  assumption  that  the  distribution  of  nucleating 
centers  throughout  the  core  volume  is  specified  by  a  Poisson  distribution.   Thus, 


>(n,V)   =  M£   e-qV  (2A) 


is  the  probability  of  exactly  n  nucleating  centers  occurring  in  a  volume,  V, 
where  q  is  the  expected  number  of  nucleating  centers  per  unit  volume.   A 
number,  k,  is  then  defined  so  that 

qV  =  k3s3  (2.5) 


1.   Haynes'  <s>  is  equivalent  to  Menyuk  and  Goodenough's  <r>.   All  of  these 
assumptions  except  (k)   and  (5)  are  explicitly  or  implicitly  contained  in  the 
equations  of  Menyuk  and  Goodenough  which  are  employed  in  the  subsequent 
development . 
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where  s  is  the  length  of  the  semi-minor  axis  of  an  ellipsoidal  domain  measured 
at  the  center.   Haynes  then  proceeds  to  compute  the  rate  at  which  core  cross- 
sectional  area  is  added  to  the  new  direction  of  magnetization.   Substituting 
the  expression  for  this  rate  into  the  equations  of  Menyuk  and  Goodenough,  there 
is  obtained  the  result 

<£  =  <cos  9>  <£>     (1  -  2e"qV)  (2.6) 

where  the  initial  condition 

§o=      -  <cos  ©>^s  (2.7) 

has  been  included.   In  the  above  expressions 
<£    is  the  magnetic  flux  and 

/£)   is  the  spontaneous  or  remanent  value  of  flux. 
Normalizing  the  flux  by  defining 


X   = 


<cos  0> 
s 


(2.8) 


and  using  relation  (2.5), 

3   3 
x  =  1  -  2e"k  <S>  (2.9) 

which  can  be  differentiated  and  combined  with  equations  (2.1)  and  (2.3)  to 
give: 

-f~  =  —^     <cos  0>  (H  -  Ho)(l  -  x)(-ln  ±-f*    )2/3.   (2.10) 

It  should  be  noted  here  that  x  lies  in  the  range  (-1,  l)  with  x  =  -1 
(as  implied  by  equation  (2.7))  as  the  initial  value. 

Equation  (2.10)  can  be  written  in  terms  of  experimentally  determinable 
constants  by  considering  the  solutions  of  (2.10)  when  H  is  independent  of  time, 
i.e.,  a  step  function.   The  family  of  curves,  dx(t)/dt,for  the  parameter,  H, 
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exhibits  a  distinct  peak  value  and  peaking  time  which  allows  k  and  <s>  to  be 
related  to  S  as  defined  in  equation  (2.2).   Haynes '  final  result  is: 

-f-  =  ±M-   (H  -  H0)(l  -  x)(-ln  4-i)2/3  .  (£.11) 

w 

This  equation  will  be  considered  in  detail  in  the  next  chapter. 

II. 6  Lindsey' s  Model 

Lindsey  (1959 )>  starting  from  the  Me nyuk- Good enough  results,  also 
developed  an  equation  for  the  function  F(<r>)(see  equation  (2.3)).  The  only 
basic  difference  between  Lindsey' s  and  Haynes'  development  is  that  Lindsey 
considers  the  flux  reversal  to  be  accomplished  by  expanding  cylindrical  rather 
than  elongated  ellipsoidal  domains.   Lindsey  finds  that  the  core  output 
voltage  for  constant  applied  driving  field  is  given  by: 

■& .  ■  h  *  ^  ■  <"■> 

whereas  the  equivalent  form  of  Haynes '  result,  as  will  subsequently  be  shown, 
is : 

k_HJt  e  .  (2.13) 


dt       3  e 
These  results  will  be  compared  in  the  next  chapter. 
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III.   FULL  SWITCHING  EQUATIONS 
III.l  Introduction 

In  order  to  provide  a  basis  for  extension  of  a  model  of  switching 
into  the  partial  switching  mode,  it  is  desirable  to  determine  which  model 
of  full  switching  most  closely  corresponds  to  reality.   This  is  done  in  the 
present  investigation  by  generalizing  the  results  of  Haynes  and  Lindsey  and 
examining  the  ensuing  family  of  equations. 

The  generalized  full -switching  equation  is  shown  to  have  certain 
mathematically-desirable  properties:   first,  the  product  of  peaking  time, 
t  ,  and  normalized  peak  amplitude,  x  ,  is  independent  of  driving  field, 
threshold  field,  and  switching  constant;  second,  the  experimentally-found 
curvature  of  the  reciprocal  switching  time  curve  can  be  duplicated  by  varying 
a  single  parameter,  n;  third,  any  usual  definitions  of  switching  time  are 
also  shown  to  be  directly  proportional  to  one  another,  in  an  easily  cal- 
culated way,  for  any  member  of  the  equation  family. 

On  the  latter  point,  it  should  be  noted  that  many  definitions  of 
switching  time  have  been  employed.  For  example,  the  time  during  which  the 
switching  waveform  amplitude  exceeds  10$  of  its  peak  amplitude  was  used  by 
Menyuk  and  Goodenough  (1955).  The  time  for  90$  of  flux  switching  has  also 
been  used,  e.g.,  by  Shevel  (1958). 

It  is  also  shown  that  within  the  family  of  equations  considered, 
the  introduction  of  radial  effects  (integration  of  the  curves  over  radius) 
does  not  lead  to  the  curvature  of  the  inverse  switching  time  curve  which 
is  experimentally  observed. 


1.   Haynes  pointed  this  out,  essentially,  for  his  equation  (n  =  3). 
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III. 2  Generalized  Full  Switching  Equations 

Haynes '  equation  is  analytically  soluble  for  H  independent  of  t,  i.e., 

for  a  step  function  of  applied  driving  field.   Under  this  condition  and  for 

H  -  H  >  0.  it  is  found  that: 

o    ' 

x  =  1  -  2  exp[-(Kt)3]  (3.1) 

x  =  6K3t2  exp[-(Kt)3]  (3.2) 

4.82  H  "  Ho 
where  K  =  —5 —  — = 

w 

x  =  — <—     where  0  =  (ft      <cos  9>,  the  observed  remanent  flux 

0  s    ^s 

s 

dx 

x  =  IT 
dt 

H  =  the  net  applied  field  intensity 

H  =  the  threshold  field 
o 

S  =  the  switching  constant,  and 

w  °  ' 

t  =  time. 
It  has  been  assumed  here  that  x  =  -1  for  t  =  0  and  that  x  has  the  range 
[-1,  l].   Note  that  H  -  H  >  0  must  hold  for  the  meaningful  application  of 
the  equations. 

We  shall  now  consider  a  family  of  equations  of  which  (3«l)  is  a 
member,  namely, 

=1-2  exp[-(Kt)n]  (3.3) 

n  >  1  being  the  region  of  interest.  Then  differentiating  (3. 3)  with  respect 
to  t, 

x   =  2nKntn'1  exp[-(Kt)n]  .  (3JO 


x 

n 
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The  time,  t      ,  at  which  the  maximum  value  of    (3»*0  occurs   is  given  by  the 
solution  of: 

x       .     0     -     SnK11  ft^C-xrfV1-1]   +   (n  -   l)tn"2  }  exp    /  -  [Kt      ]n  }        . 

n  L  np       np     v       np  J    r  1     np  J 

-nKV  +  (n  -  1)  =  0 
np 

i^i) 1/n 

%  ■   ^  V         ■  (3.5) 

Substituting  (3.5)  into  (3«*0i  and  setting  x  (t   )  =  x  , 
DV  v-*  "  nv  np       np' 

/n  -   1\  (n'l)/n 

i         =     2nKn     N     n      )    exp    /-     S^l\ 

np  jjn-1  I  n      J 

, .     (n-l)/n  

i  =     2n     5_Z_1\  e  n       K  .  (3.6) 

np  \_     n      J 

Multiplying    (3.5)  by  (3-6),   we  obtain 

n-1 

t     ±  =     2n  (*LZ-±)     e   ' 

np  np 


*m 


n-l 


t     x  =     2(n  -   1)   e  n  .  (3.7) 

np  np  v  v 

Hence,  we  have  the  conclusion  that  x  t   is  dependent  only  on  n, 

np  np  ' 

being  independent  of  K,  regardless  of  the  possible  dependency  of  K  on  n, 

Also,  the  value  x(t  )  =  x   is  independent  of  K,  which  can  easily 

be  seen  by  substituting  (3. 5)  into  (3«3)>  whereupon: 

n-l 
xp  =  1  -  2e    n    .  (3.8) 
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It  can  also  be  shown  that  any  given  fraction,  a,  of  the  peak  time 
rate  of  flux  change,  x  ,  occurs  at  a  time  which  is  a  constant  multiple,  P, 
of  the  peaking  time,  independent  of  the  value  of  K.  Thus,  using  (3.^)> 

X  =      OK  =       2rfSn-VClEfc]n 

n      np 
and  setting  t  =  pt    where  0  <  a  <   1  and  p  >  0 

n   ^Qn-1  n-1  -[Kpt   ]n 
ax         =     2nKT  6   t    e   ^  np 
np  np 

Substituting  (3*5)  and  (3.6), 

,    lx  (n-l)/n    -£=i  n   ,       (n-l)/n   -Pn(^) 

a  .  2n  (H^-I)         e    n  R  =  2npn-lK  ^_lj         e   U  j 

n-1    ¥  ^^ 
a  =  P     e  (3.9) 

which  establishes  the  preceding  statement.   Two  solutions,  f3,  and  |3  ,  one 
before  and  one  after  peak  amplitude,  exist  for  each  0  <  a  <  1  for  a  given 
n. 

Any  given  value  of  x  can  also  be  related  to  a  constant  (relative 
to  K)  by  substituting  t  =  pt   into  (3«3)>  from  which: 

n-1   n 

n~  P 

x  =  1  -  2e  .  (3-10) 

Therefore,  for  any  member  of  the  family  of  equations,  i.e.,  for 
any  n  >  1,  the  relation  of  x  or  of  any  fraction  of  x  to  the  peaking  time 
is  constant.   Hence,  all  usual  definitions  of  switching  time,  whether  given 
by  the  time  to  some  value  of  x  or  some  fraction  of  x  ,  are  conveniently 
interrelated  within  this  family.   Similary,  the  ratio  of  switching  time  to 
peaking  time  is  a  constant  for  any  n.   The  meaning  of  non- integer  n  is 
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problematical  and  will  be  discussed  in  a  more  appropriate  context 
(Chapter  V). 

III. 3  Fitting  Experimental  Data 

Table  1  and  Fig.  3.1  present  a  typical  group  of  experimental  values 

of  x  t   (the  peak  amplitude -time  product),  which  are  observed  for  full 

switching,  i.e.,  for  a  saturated  state  as  the  initial  condition  of  the  core. 

The  value  of  x  t   is,  clearly,  a  monotonically  increasing  function  of  F, 

the  applied  mmf .   Since  the  experimental  value  of  x  t   is  not  constant,  the 

model  of  full  switching,  which  predicts  constant  x  t  ,  is  limited  in 

'  p  p 

accuracy,  for  any  single  n  in  equation  (3«3)>  at  the  outset.   It  should  be 
noted  that  since  equation  (3-7)  is  independent  of  K,  it  is  impossible  to 
contrive  a  function,   K(H,  H  ,  S  )  which  conforms  to  experiment  over  a 
range  of  H  (with  n  constant  and  H,  H  ,  and  S  independent  of  time).  There- 
fore, both  Haynes'  and  Lindsey's  models  (which  are  members  of  the  class  of 
equations  being  studied)  are  faulty  to  some  degree. 

Equation  (3*7)  is  plotted  as  Fig.  3.2.  The  values  of  n  which 
correspond  most  closely  to  the  various  experimental  values  of  x  t  may  be 
conveniently  found  from  Fig.  3 .2.  The  values  so  obtained  vary  from  about 
1.8  to  2.2  for  the  range  of  applied  drives  listed  in  Table  1,  which  is 
the  range  of  interest  for  later  work.   Hence,  the  vicinity  of  Lindsey's 
equation  (n  =  2)  corresponds,  on  this  basis,  to  the  more  nearly  correct 
model . 

III. k     Curvature  of  the  1/T    versus  F  Curve 
[ s 

One  of  the  basic  relations  of  the  Menyuk-Goodenough  (1955)  model 

is,  for  step  functions  of  H  and  assuming  initial  saturation, 


1.   One  case  does  not  conform,  this  being  attributed  to  "experimental  error" , 
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Table  1 

F 
(ma.-t. ) 

0 

(mv. ) 

XP   '"  Js 
(us.)" 

t 
P 

(|1S.) 

x  t 
P  P 

450 

71.5 

3-03 

0.29 

O.878 

500 

9^-5 

4.00 

0.24 

O.96 

550 

123 

5.22 

0.19 

0.992 

600 

150 

6.36 

0.17 

1.082 

650 

178 

7.5^ 

0.l4 

1.057 

700 

208 

8.82 

0.13 

1.148 

750 

240 

10.17 

0.118 

1.198 

800 

27^ 

11.60 

0.105 

1.218 

850 

308 

13.03 

0.095 

1.238 

0        =     23.6    •    10 "^  webers 
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(H  -  H  )  ■£-     =     S 
v     o   T  w 

&  s 

where   £*    is  the  switching  time  and  is  defined  as  the  interval  between 

5 

the  points  of  10$  peak  amplitude  (of  the  irreversible  component). 
This  relation  was  inferred  from  experiment  and  is  often  presented  as  a 
linear  relation  between  l/f  t   and  F,  the  applied  mmf .   It  should  also  be 

5 

noted  that  if  the  abscissa  of  such  a  curve  is  labelled  in  units  of  field 
intensity,  this  is  always  obtained  by  some  relation  of  the  form 

F 

H   =  T 

xm 
JL     being  the  mean  path  length  through  the  magnetic  specimen.   This  curve  is 
usually  called  the  "inverse  switching  time"  curve. 

It  was  found  experimentally  by  Gyorgy  (1957)  and  later  verified  by 
many  investigators  using  fast  rise-time  drive  pulses,  that  the  inverse 
switching  time  curve  is  concave  upward.   Gyorgy  used  a  different  definition 
of  switching  time  (viz.,  the  time  from  10$  to  90$  of  flux  change)  than  most 
earlier  investigators,  but  there  is  no  reason  to  doubt  that  his  definition 
differs  by  other  than  a  constant  from  the  earlier  definition. 

Fig.  3«3  is  a  graph  of  v  (the  peak  amplitude  of  the  irreversible 

flux  change)  plotted  against  drive  mmf,  F.   It  is  also  concave  upward,  which 

reflects  the  fact  that  since  \j%       (or  l/t  )  is  increasing  faster  than 

linearly  and  the  area  under  the  switching  waveforms  is  constant,  the  peak 

amplitude  is  likely  to  rise  at  a  faster-than-linear  rate.   This  effect  is 

also  seen  in  the  increasing  values  of  x  t   in  Table  1. 

P  P 

It  is  therefore  not  surprising  that  neither  Haynes '  nor  Lindsey's 
equation  exactly  fits  experimental  observations,  since  both  are  based  on  the 
linear  reciprocal  switching  time  relation  through  the  Menyuk-Goodenough 
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model.  On  the  other  hand,  either  of  the  equations  is  of  sufficient  accuracy 
for  a  number  of  applications. 

A  further  point  should  be  made  concerning  the  deduction  of  physical 
facts  from  the  model(s).  The  results  of  the  preceding  section  show  that  no 
function  of  the  parameters  H,  H  ,  and  S   (independent  of  t )  could  cause  the 
observed  non-zero  curvature  of  the  inverse  switching  time  curve.   It  would 
be  tempting,  but  fallacious,  to  conclude  that  since,  for  example,  H  -  H 
is  assumed  proportional  to  the  domain  wall  velocity  for  the  model,  then 
(physical)  variations  of  domain  wall  velocity  cannot  be  the  cause  of  the 
observed  curvature. 

The  accuracy  of  the  present  model  is  too  poorly  established  to 
allow  any  degree  of  assurance  in  predicting  physical  phenomena. 

III. 5  Radial  Effects 

During  flux  reversal  in  square-loop  cores,  the  average  permeability 
varies  widely  with  the  transient  state  of  magnetization.   For  a  toroid  of 
uniform  high-permeability  magnetic  material,  the  magnetizing  flux  is  nearly 
uniformly  distributed  over  any  cross- section  in  the  material  perpendicular 
to  the  flux  lines.   That  is,  since   B  =  u  H  +  M  and,  in  this  case, 
| M |  »  |u  H|  ,   the  variation  of  H  with  path  length  (in  agreement  with 
Ampere's  law,   F  =  9   H  ■  dl)  is  masked  by  M. 

If  the  permeability  is  low,  however,  |m]  is  of  the  order  of  |h| 
so  that  the  flux  is  more  non-unif ormly  distributed  over  the  core  cross- 
section.   Since  the  field  then  depends  on  radius,  this  state  of  affairs  is 
called  a  radial  effect. 

During  switching  of  a  square-loop  core,  it  is  clearly  quite 
difficult  to  determine  the  inaccuracy  involved  in  assuming  a  uniform  flux 
distribution. 
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Although  appreciable  radial  effects  do  occur  experimentally,  it 
is  expected  that  this  variation  of  flux  density  over  core  radius  can  be 
simulated  by  an  expeditious  choice  of  the  statistical  distribution  of  H 
(or  a  related  quantity  such  as  domain  vail  velocity). 

Two  groups  of  investigators  reported,  nearly  simultaneously, 
results  dealing  with  radial  effects.   Tancrell  and  McMahon  (1959)  assumed 
a  simple  waveform  for  the  switching  of  a  thin  annulus  of  ferrite  material, 
this  being  of  the  form  of  a  sinusoid.   Integration  was  performed  over 
radius  for  various  values  of  parameters,  which  yields  curves  closely  re- 
sembling a  full  switching  curve.   Bardezh  and  Kobelef  (1959)  report  very 
similar  results.   Tancrell  and  McMahon  also  confirmed  their  results  by 
drilling  small  (.005")  sensing  holes  at  various  radii  in  a  thick-walled 
core  and  observing  the  waveforms  of  the  various  layers.   The  experimental 
and  theoretical  switching  times  were  said  to  agree  with  calculated  values 
to  within  10%. 

These  results  strongly  support  the  case  for  the  existence  of 
radial  effects. 

If  radial  effects  must  be  treated,  then  a  statistical  approach 
would  seem  to  offer  an  attractive  alternative  to  the  lengthy  integration 
otherwise  required.   Using  the  statistical  approach,  there  would  remain 
some  question  of  the  accuracy  with  which  true  radial  effects  could  be 
duplicated. 


III. 6  Radial  Models  for  Haynes '  and  Lindsey's  Equati 


ons 


If  Haynes'  equation  is  considered  as  the  equation  for  switching 
of  a  differential  thickness  of  core  material,  the  switching  waveform  for 
the  entire  core  is: 


*L  M     £n£\3       o  > 
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0  =  v7^''6!— J  *  y  0-H°)  exp  -hr(/-H° 

where  Jt  ,  X,   are  "the  inner  and  outer  path  length  respectively. 

The  above  integration  yields  data  typically  as  summarized  in  Fig.  3«^«  The 

result  of  using  Haynes '  equation  for  the  total  core  assumed  centered  at 

mean  radius  is  shown  as  a  reference.   In  all  cases,  H  and  S  are  assumed 

o      w 

to  be  constants. 

It  can  be  seen  that : 

(1)  Integration  over  radius  produces  no  curvature  of  the  inverse 

switching  curve  or  of  the  v  versus  F  curve  (as  long  as  F 

P 

is  large  enough  to  switch  the  outer  radius). 

(2)  An  effect  of  radial  integration  is  to  decrease  the  peaking 

time,  t  ,  and  the  peak  amplitude,  v  ,  relative  to  the 

result  of  using  the  same  H  and  S  in  Haynes'  equation 

o      w 

directly. 

(3)  After  integration,  the  F  axis  intercepts  of  the  l/t  and ' 

the  v  curves  are  no  longer  identical. 
P 

These  results  have  proved  to  be  typical  for  integration  over 
equations  of  the  family  defined  by  (3^)-  In  particular,  this  is  true 
for  Lindsey's  equation. 

Result  (l)  above  therefore  tends  to  support  the  argument  that  the 
curvature  of  the  inverse-switching-time  curve  is  not  due  to  radial  effects 
and  therefore  implies  a  material  effect  as  the  cause.   The  nature  and 
extent  of  the  various  assumptions  which  have  been  made  in  the  course  of 
the  development  do  not  permit  any  broad,  definite  conclusions. 


d* 


A  Haynes '  Eqn.  unaltered  (mean  path  length  =  3  nra) 

D  Haynes'  Eqn.  integrated  over  path  lengths 
(2.238  -  3.77  mm) 
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Result  (2)  shows  its  effect  in  a  larger  and  longer  tail  of  the 

switching  waveform.  That  is,  since  t  is  smaller  and  v  is  smaller,  the 

P  P 

constancy  of  area  under  the  switching  curve  requires  a  slower  decay  of  the 

curve . 

The  value  of  x  t  resulting  from  integration  can  be  matched  within 
p  p        °         ° 

the  family  of  curves  defined  by  (3*^)  by  choosing  a  smaller  value  of  n  to 
represent  the  total  core  than  was  used  for  material  alone. 

Result  (3)>  however,  shows  that  the  integration  of  a  member  of 
(3.^)  over  radius  produces  a  curve  which  is  not  a  member  of  (3-^)  since 
the  members  of  (3.^)  have  identical  F-axis  intercepts  for  the  l/t  and  v 
curves.   Therefore,  it  is  not  possible  to  exactly  account  for  radial  effects 
by  adjustment  of  n.   It  is  not  to  be  expected  that  equation  (3^)  can  yield 
an  arbitrarily  accurate  switching  waveform  for  a  core  of  any  desired 
radial  ratio. 
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IV.   THE  MAXPEAK  EFFECT 
IV. 1  Introduction 

It  was  demonstrated  in  the  previous  chapter  that  the  full  switching 
models  are,  to  some  degree,  inexact.  The  attempt  to  locate  the  origin  of 
some  of  the  inaccuracy  led  to  the  investigation  of  partial  switching,  that 
is,  of  flux  changes  initiated  from  a  non- saturated  state.   It  was  expected 
that  the  additional  information  available  in  partial  switching  experiments 
would  provide  more  penetrating  information  concerning  domain  activities  than 
would  full  switching  experiments.  Thus,  the  attempt  to  construct  a.  model  of 
a  partial  switching  effect  could  be  expected  to  shed  some  light  on  the  inac- 
curacies of  full  switching  models. 

In  order  to  be  useful,  the  selected  partial  switching  effect  should 
be  quite  marked  and  experimentally  reproducible.  The  "maxpeak  effect"  was 
found  to  meet  these  requirements  very  well. 

The  construction  of  a.  model  of  the  "maxpeak  effect"  is  the  prin- 
cipal objective  of  the  following  investigation.   In  developing  this  model, 
we  shall  endeavor  to  retain  a  correspondence  between  the  mathema.tics  and  the 
physical  concepts  involved,  rather  than  proceed  along  strictly  mathematical 
lines.   It  will  be  shown  that  a  potentially  useful  model  of  the  maxpeak 
effect  can  be  developed  by  logical  extension  of  the  concepts  of  domain 
activity  described  by  Manyuk  and  Goodenough  (1955)  and,  in  more  detail  by 
Haynes  (1958)  and  Lindsey  (1959);  particularly  the  latter. 


1.  The  usual  meaning  of  "partial  switching"  is  somewha.t  broa.der  than  the 
present  usage.   Usually  any  flux  change  substantially  less  than  the  maximum 
possible  (twice  remanent  flux)  is  called  partial  switching,  regardless  of 
the  initial  state. 


3^ 


Evidence  is  presented  which  supports  the  idea,  that  a  considerable 
distribution  of  domain  wall  velocities  (for  given  driving  field)  may  be  in- 
volved in  the  flux  reversal  process.   It  is  shown  that  the  introduction  of 
such  a.  velocity  distribution  is  harmonious  with  other  relevant  data,. 

Several  negative  results  of  the  search  for  a  model  are  also 
given. 

In  the  following  description  of  experiments,  the  reproducible 
saturated  state  of  the  core  is  called  the  DOWN  state.  When  the  core  is 
driven  toward  the  other  saturated  state,  it  is  said  to  be  switched  or 
driven  UP  to  the  UP  state . 

IV. 2  Definition  of  the  Maxpeak  Effect 

Consider  a  magnetic  core  to  be  driven  by  a  series  of  pulses  a.s 
indicated  in  Fig.  k-.l   with  the  core  set  to  the  DOWN  state  initially. 


450  to  1100  ma.-t. 
variable  amplitude 

UP   i , 

400  mus. 

typical 


UP 


DOWN 


UP 

200  mjis., 
typical 

(fixed  for 

each 
experiment) 


t 


RESET 
(DOWN) 


DN 


i+50  to 
700  ma.-t 


2.5  to  3  a.-t 
15  us. 


Figure  k.l 


As  the  value  of  F   is  increased,  the  core  is  driven  successively 
closer  to  the  saturated  UP  state  and  the  sequence  of  successive  switching 
waveforms  due  to  the  DOWN  drive  is  shown  in  Fig.  k.2   (see  also  Section  IV. k) 
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It  can  be  seen  that  the  pea.k  amplitudes  of  the  curves  go  through  a  maximum, 
in  this  case,  for  the  amplitude  F   =  800  ma.  This  peak  amplitude  is  termed 
the  "maxpeak"  while  the  phenomenon  is  generically  named  "the  maxpeak  effect". 
This  effect  is,  visibly,  a  substantial  one,  the  maxpeak  amplitude  being  some 
30  -  90$  greater  than  the  full-switching  peak  amplitude  (for  the  same  F  ). 

For  the  present,  we  shall  only  consider  the  maxpeak  effect  at 
higher  values  of  F   (corresponding  approximately  to  2H  £=J  700  ma.  -t). 

In  all  subsequent  work,  rectangular  current  pulses  having  rise 
times  of  1  mu.sec.  (hence  only  0.5$  of  the  average  switching  time)  are  used. 
The  RESET  drive,  for  which  only  peak  amplitude  and  duration  are  critical, 
has  a  somewhat  slower  risetime  (see  Appendix  i). 

IV. 3  Previous  Investigations 

Smith  (i960)  is  the  only  investigator  known  to  have  published  any 
data  concerning  the  maxpeak  effect.   His  experiments  were  performed  on  a 
specially  prepared  thin  -  wall  core  having  an  OD/lD  ratio  of  1.2.   The 
material  of  the  core  is  probably  very  similar  to  that  of  the  core  used  in 
the  present  study,  a  MgMn  ferrite  (see  Appendix  i).   The  conclusions  which 
Smith  reaches  concerning  this  effect  are: 

(1)  the  difference,  (maxpeak  amplitude) -(full  switching  peak  amplitude) 
=  a.  constant,  viz.,  60  mv.  over  a  wide  range  of  drives. 

(2)  the  maxpeak  curve  contains  90$  of  the  total  flux  of  the  core 
i.e.,  maxpeak  occurs  for  the  condition  that  90$  of  the  total 
flux  wa.s  switched  during  the  UP  pulse  in  all  ca.ses. 

His  investigation,  not  being  concerned  with  models,  proceeds  only 
to  the  point  of  stating  the  above  observations  and  some  general  comments  on 
the  possible  origin  of  this  effect. 
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IV .  ^4-  Experimental  Results 

Fig.  U.3  presents  curves  of  the  pea.k  amplitudes  of  switching 

curves  under  various  conditions  of  F   and  T   for  a,  given  F   .   It  is 

immediately  seen  that  the  difference  v      ,  -  v _  , .    . ,  .  .      ,  is 

maxpeak    full  switching  peak 

remarkably  close  to  the  value  given  by  Smith.   In  the  present  case,  the 
values  are  69,  lh,    and  58  mv.  (a,verage  67  mv.)  respectively  as  compared  to 
Smith's  value  of  60  mv.,  in  spite  of  the  difference  in  OD/lD  ratio  of  the 
cores  used  in  the  two  investigations  (Smith  1.2;  present,  1.68). 
We  now  introduce  the  symbolism  defined  by  Fig.  h.k. 


or  v 


maxpeak  switching  curve, 
v 


ull  switching  curve, 
*fs  or  vf s 


Figure  h.k 


From  the  experimental  data,  we  deduce  the  following  statements 
concerning  the  maxpeak  switching  curve. 
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(l)         fxx  at    *    0.83    /   Xj 


CO 


o 


^fsdt 


(2)  t-^E     £f     1.3  near  Fn.T  =  2H     =700  ma. 

x„  DN  o 


fsp 


in  present  case. 


(3)       t    =  0.88  t_  fh   .. 

v  '  xp  fsp  (^.l) 


It  should  be  noted  that  for  (l)  above,  Smith  would  ha.ve  0.9  rather 

than  O.83.   Since  Smith  concludes  that  x   -  x_   =  constant,  there  is  no 

xp    fsp  ' 

exact  equivalent  to  (2)  above.  Also,  no  relation  of  the  variables  of  (3) 

above  is  given  by  him. 

In  a  later  section,  evidence  will  be  presented  which  demonstrates 

that  neither  relation  (2)  or  (3)  holds  over  extensive  ranges  of  F   and  that 

Smith's  relation  x   -  x„   =  constant  is  likewise  inadequate. 

xp    fsp 

IV. 5  The  Formulation  of  a  Model 

The  object  of  the  following  investigation  is  to  discover  a.  model 
of  switching  which  describes  the  maxpeak  effect.  The  principal  constraint 
in  this  inquiry  is  that  the  model  be  based  on  switching  processes  which  are 
either  known  or  are  reasonable  on  the  basis  of  present  evidence.   Stated 
negatively,  purely  mathematical  models,  having  no  known  relation  to  physical 
factors,  are  specifically  disallowed. 

It  is,  of  course,  desirable  that  any  such  model  be  as  simple  as 
possible  within  certain  accuracy  limitations.  A  case  in  point  is  the 
lingering  question  of  radial  effects.   If  they  must  be  introduced  exactly 
the  mathematical  complexity  is  greater  due  to  the  lengthy  integration 
which  ensues.   It  is  desirable,  on  the  other  hand,  to  account  for  radial 
effects  statistically.  Both  the  exact  and  statistical  approaches  are 
subject  to  consideration,  however. 


IV. 6  The  Maxpea.k  Effect  on  the  Basis  of  Ra.dia.l  Effects  Only 

Prior  to  the  publication  of  Smith's  work,  the  possibility  that 
the  maxpeak  effect  is  due  solely  to  radial  effects  was  considered.  For 
this  investigation,  switching  curves  of  the  family  of  equations  (3.3)  and 
(3.^0  were  used.   It  is  assumed  that  the  switching  curve  followed  by  a 
given  annulus  of  the  core  is  the  final  portion  of  the  full -switching  curve 
which  would  have  been  followed  if  switching  had  commenced  at  saturation. 
The  following  equations  apply. 
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of  UP  drive  -(KT    )n 


X  l-2eUU  t       =     T 

u  u  u 


.      -[K.(t.  +  T  )]n 
_     „n  /.         /r-\n-l  dx    d  ' 

xd     =      -2n  K^  (t     +T) 


-[Kd(t   +  r)]n 

During  DOWN  Drive    xn   =   -1  +  2e 
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.      rr*        1  /  ,       U\l/n        t,  >  0 

where  Z  =     =  ( -In  —^ )  ;  d  - 

Kd       2 


Subscript  u  corresponds  to  UP,  d  to  DOWN.   Capitals (X,  X,  T) 
imply  final  values. 

These  equations  are  then  integrated  (see  Appendix  II)  over  values 
of  K  and  K,  corresponding  to  the  various  rings  into  which  the  core  is 
divided.  By  adjusting  H  ,  S  ,  a.nd  n,  it  is  possible  to  obtain  results 
similar  to  those  of  Fig.  4.3. 


in 

Fig.  4.5  is  a  case  in  which  the  curve   for  n  =  3  provides  a 

remarkably  good  fit  to  the  experimental  data.   In  all  cases,  the  values 

of  H  and  S  are  adjusted  such  that 
o      w 

(1)  the  absolute  maxpeak  amplitude  is  correct  and 

(2)  the  maxpeak  occurs  at  the  correct  value  of  F , 
the  driving  m.m.f.  in  the  DOWN  direction. 

The  t  curves  for  n  =  3  and  n  =  k   are  quite  tolerably  near  the 
experimental  curve  while  for  n  =  2,    the  agreement  is  poor.  Also  for  n  =  3, 
the  absolute  full-switching  peak  amplitude  is  in  good  agreement.   However, 
when  this  model  is  applied  to  the  case  of  maxpeaking  arising  from  other 
values  of  T  ,    it  is  found  that  the  required  value  of  H  becomes  physically 
unreasonable.  This  fact  alone  leads  to  the  conclusion  that  the  maxpeak 
effect  is  probably  not  physically  ascribable  to  radial  effects  only. 

Further  confirmation  of  this  conclusion  is  given  in  Section  IV. 9 • 

IV. 7  Basic  Model  of  Maxpeak  Effect 

From  the  similarity  of  Smith's  results  with  the  present  results, 
as  indicated  in  Section  IV. k,    it  is  logical  to  a.ssume  that  the  maxpeak 
effect  depends  principally  upon  microscopic  magnetic  properties  rather 
than  macroscopic  ones.  An  extension  of  Lindsey's  model  to  explain  this 
effect  will  now  be  considered.   (Lindsey's  model  is  used  because  the  results 
of  Section  III. 3  indicate  n  =  2  agrees  best  with  full -switching  data.) 

The  equation  of  the  full -switching  curve  for  n  =  2  may  be 
written: 

x  =  4K2t  e_ICt  (4.3) 
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^3 
Examination  of  the  derivation  of  this  equation  shows  the  following  signif- 
icance of  the  various  terms: 

K  is  the  domain  wall  velocity  in  sec.   ;  that  is,  the  velocity 
normalized  to  length  units  of  core  size;  or,  if  preferred, 
l/K  is  time  required  for  the  fraction,  1-e  ,  of  the  flux  to 
be  reversed. 
k&t     is  the  time  rate  of  transfer  of  cross -sectional  core  area 

( /v/  flux)  to  the  new  direction  if  no  domain  wall  collisions 

occur. 

.j{2t2 

e      is  the  probability  of  a  domain  wall  having  not  yet  been 

annihilated  at  time  t;  or,  the  fraction  of  nucleating  centers 
whose  walls  are  still  "active"  at  time  t. 

Consider  a  cross -section  of  polycrystalline  material  having  been 
switched  such  that  the  domain  walls  are  as  indicated  in  Fig.  k,6a..      Some 
domains  have  stopped  growing  (according  to  either  Lindsey's  or  Haynes ' 
model)  at  small  radius,  others  at  larger  radius  (type  a)  and  still  others 
(type  p)  would  continue  to  grow  if  the  field  were  not  removed  at  time  T   . 
Upon  application  of  a  sufficiently  large  drive  field  in  the  direction  H-^, 
it  is  assumed  that 

(a)  inactive  domains  (type  a  domains)  simply  renucleate 
from  their  original  nucleating  centers. 

(b)  active  domains  (type  3  domains)  both  renucleate  from 
their  original  centers  and  the  existing  domain  wall 

is  also  pushed  inward,  preserving  its  circular  geometry. 
Domains  in  this  class,  when  subjected  to  the  DOWN  pulse 
will  deliver  flux  in  a  way  shown  by  the  following  analysis. 


kk 


Referring  to  Fig.  k.6,    if  R  is  the  radius  of  the  domain  at  time 


•  "up  or  V 

+  HDN°r  * 


DN 


(a) 


00 


Figure  k.6 

T  ,  then  at  some  time  during  the  DOWN  pulse,  the  radius  of  the  outer  wall 
will  be  (R-r)  while  the  newly  renucleated  inner  wall  is  at  radius  r  since 
the  velocity  of  inner  and  outer  walls  is  equal.  The  total  length,  X.  ,  of 
domain  wall  at  this  time  is: 


$  =  2irr    +  2fT(R-r)  =  2rrK 


dA 


CLA.      0  CLA. 

Since  —  =  vX  ,   where  v  is  the  wall  velocity,  it  can  be  seen  that  —  and, 

therefore,  -r£-  is  constant,  for  this  type  of  domain  activity  until  the  inner 

■p 
and  outer  walls  collide.  This  latter  condition  corresponds  to  r'  =  —  which, 

since  R  ^  KT  while  r'  =  K,t'   can  be  written, 
u  u  d  d'  ' 


Vi 


K  T 
u  u 


[h.k) 


t*  being  the  time  of  annihilation  of  the  two  counter -moving  (3-domain  walls 


Meanwhile,  the  walls  renucleated  from  a-type  centers  in  Fig.  4.6a,  simply 
retrace  the  flux  reversal  procedure  of  the  previous  UP  pulse. 

IV. 8  Analysis  of  the  Basic  Model 

The  model  proposed  for  the  maxpeak  effect  can  now  be  conveniently 
examined  graphically.   (The  exact  description  is  given  by  equations  (4.5)). 

The  details  of  this  graphical  method  will  be  demonstrated  by  the  use  of 

2 
Fig.  4.7.   Since  the  quantity  K  /K,   is  only  a  scaling  factor  of  the  amplitude 

and  time,  setting  K  /K,  =  1  greatly  simplifies  the  presentation. 

In  Fig.  4.7,  the  two  components  of  flux  which  have  been  reversed 
at  time  T  are  indicated  as  the  shaded  areas  in  the  upper  diagram.   The 
flux  transfer  during  the  DOWN  drive  is,  of  course,  the  chief  concern.  The 
amplitude,  x,,  is  actually  of  opposite  sign  to  x  but  is  shown  the  same  for 
convenience.   During  the  DOWN  drive,  the  flux  change  is  hypothesized  to  con- 
sist of  two  parts,  the  a   and  the  (3  domain  contributions.  The  former  contri- 
bution has  precisely  the  same  waveform  during  the  UP  and  the  DOWN  drive 

(remembering  the  assumption  K  =  K,). 

u    d 

The  contribution  of  the  8  domains  is,  itself,  composed  of  two 
parts.  The  inner  walls  (renucleated  walls)  contribute  a,  component  of  flux 
exactly  equal  to  the  contribution  of  8  domains  from  time  t  =  0  to  t  =  T  /2 
(see  below)  during  UP  drive.  The  outer  6  domain  walls  (those  already  nucle- 
ated and  not  annihilated  during  the  UP  drive)  contribute  the  reaminder  of 
flux  delivered  by  8  domains .  The  total  time  rate  of  flux  change  during 
DOWN  drive  is  the  sum  of  the  a   and  8  contributions  as  indicated  in  Fig.  4.7. 

At  time  T  /2  =  t',  the  inner  and  outer  8  domain  walls  collide 
u'     d' 

since,  as  shown  in  the  previous  section,  t'  =  K  T  /2K„  which  is  just  T  /2 

d    u  u'   d  u' 

when  K  =  K,. 

u    d 


2.   Subscript  "u"  refers  to  UP  drive,  "d"  to  DOWN  drive. 
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*UP  '  SN  °r  FUP  =  FDN'   n  =  2 


Full  switching  curve 
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Figure  4.7  Graphical  Construction  of  M&xpeak 
Effect  Curves  for  Basic  Model 


^7 
The  preceding  description  corresponds  to  some  single  set  of  F 

(or  K  ),  F,  (or  K, ),  and  T  .   In  order  to  obtain  the  family  of  curves  for 
v    u    d      d       u 

the  maxpeak  effect  similar  to  Fig.  h.2,    the  value  of  F  (or  K  )  is  usually 
varied  while  T  and  F  are  held  constant.   For  the  purposes  of  graphical 
analysis,  it  is  much  simpler  to  vary  T  than  F  .   It  is  possible  to  inter- 
change these  two  variables  because,  for  a  given  n,  Kt  always  occurs  as  a 
product  in  the  switching  equation  (Eq.  3.*0  except  for  a  multiplier  of  K. 
This  multiplier,  therefore,  simply  scales  the  absolute  magnitudes  of  the 
amplitude  and  time  axes .   It  has  been  tacitly  assumed  that  n  is  the  same 

for  F  and  F , . 
u      d 

The  maxpeak  effect  curves  corresponding  to  the  "basic"  model  of 
the  maxpeak  effect,  constructed  by  the  method  indicated,  are  presented  in 
Fig.  ^.8.  The  resultant  curves,  numbered  1  to  3j  are  each  constructed  as 
explained  in  connection  with  Fig.  k.J.      It  can  be  seen  that  the  maxpeak 
curve  is  number  2  whose  peak  amplitude  is  1.33  xf 

Several  differences  between  the  maxpeak  curve  deduced  from  the 
model  and  that  found  experimentally  are  apparent.  The  model  predicts 
maxpeak: 

(1)  at  a  value  of  flux  which  is  too  small  (  m  63$  rather 
than  83$  of  the  total  flux), 

(2)  at  a  time  which  is  too  early  relative  to  t„ 

fsp 

Also,  the  wall  collapse  of  B  domains  shows  a,  deterministic  fall  of  voltage 
which  is  not  observed  experimentally.   In  this  connection,  it  is  also  noted 
that  Smith's  curves  (for  a.  thin  core)  show  voltage  falls  no  more  precipitous 
than  for  a  thick  core,  suggesting  that  smoothing  due  to  radial  effects  is 
not  the  principal  factor  in  removing  the  sharp  drop  predicted  by  the  basic 
model. 
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^9 
On  further  examination  it  is  found  that  only  « k$   of  the  total 
flux  would  have  to  be  shifted  under  the  curve  given  by  the  model  in  order 
to  match  the  experimental  maxpeak  curve  (using  the  model-derived  curve 
which  includes  the  same  total  area,  as  the  experimental  curve). 

In  summary,  the  basic  maxpeak  effect  model,  which  is  a.  logical 
extension  of  Lindsey's  model,  has  several  shortcomings  when  compared  against 
the  absolute  standard  of  experimental  data.   Positive  results  are  also  present, 
the  most  outstanding  being  the  maxpeak  amplitude  ratio  x  /x„   =  1.33  ac- 
cording to  the  model  as  compared  to  1.30  to  I.36  for  various  experimental 
cases  (Fig.  ^.3)  for  relatively  large  F™. 

The  equations  for  the  basic  maxpeak  model  are  as  follows: 
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where 


-  K2^ 


x,         =     4ICt,     e 
dfs  dd 


*B»>  =  *a  i\  -  °>  =  (r)  K  =  (r)  ^  e "  -i 


K  T 
T'  u  u 

d  2  K, 

d 


and 


K     =     1.606 


fH  -  m 
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IV. 9  Modification  of  the  Basic  Model  by  Wall  Velocity  Distribution 

In  order  to  obtain  a  more  accurate  model  for  the  maxpeak  effect, 
various  modifications  of  the  basic  model  will  be  considered. 

The  deterministic  fall  of  voltage  which  the  basic  model  predicts 
is  of  immediate  concern.   Studies  (e.g.,  LFE  i960))  of  thin  metallic  films 
have  shown  that  variations  in  the  velocities  of  domain  walls  occur  both  for 
individual  walls  as  functions  of  position  and  for  all  walls  collectively. 
Such  variations  may  be  due  to  composition  gradients  of  various  kinds. 
There  is  insufficient  data  available,  however,  to  adequately  assess  the 
magnitude  of  this  effect. 

Calculations  have  been  made  (similar  to  Section  IV. 8)  in  order  to 

examine  the  effect  of  velocity  variation  on  the  maxpeak  effect  curves.   It 

is  supposed  that  the  time  of  collapse  of  the  B  domains  is  normally  distrib- 

K  T  K  T 

uted  about  time      and  with  standard  deviation   #"=  0.2  (     ) .  The 
c  K-,  c.    K-, 

d  d 

resulting  set  of  curves  is  given  in  Fig.  k.9. 
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The  corresponding  velocity  distribution  should,  for  consistency, 
also  be  introduced  into  the  full -switching  equation.  The  resulting  full- 
switching  curve  is  shown  in  Fig.  Ij-.IO.   It  is  evident  that  the  error 
thereby  introduced  is  negligible. 

A  measure  of  the  relative  accuracy  of  the  maxpea.k  model  is 

afforded  by  the  average  absolute  value  of  the  error  as  a  percentage  of 

3 
peak  amplitude.   Thus, 


A  .  i   Z 
n  i-1 


w  -  w 


w 


where  vm(t.)  is  the  theoretical  voltage  at  time  t.  and  v„,r(t.)  is  the 
T  1  i      EX  i 

experimental  voltage  at  time  t . . 

For  the  data  of  Fig.  h.ty,   A  =  6.5$  is  found  for  a,  typical  maxpeak 
curve.   By  comparison,  a  full -switching  curve  from  the  family  of  equation 
(3»*0  compared  to  a  fitted  experimental  full -switching  curve  corresponds 
toA  =  2$. 

Therefore,  the  maxpeak  effect  model  with  velocity  smoothing  is 
reasonably  accurate.   It  should  be  pointed  out,  however,  that  other  highly 
sensitive  measures  of  accuracy  of  the  model  show  that  some  problems  remain. 
Namely,  the  model  predicts  the  maxpeak  at  a  flux  level  which  is  too  small 
(approximately  70$  rather  than  83$).   Also,  t  /t    is  too  small  (0.55 
as  compared  to  0.88). 


3.  This  measure  tends  to  emphasize  voltage  errors  in  accordance  with  the 
interests  of  engineering  applications. 
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The  equation  for  the  maxpeak  model  with  velocity  distribution 
is  as  follows: 


t 


(y  -  Tp2 


*d     -     ^d     e"(V/      -     &^td  +  {ii       (.1./^      e     ^         J 

d  dd  \V     Tu     d       Ku     U     1         -i-***  J 


for  tJ     >     0 
d     — 


All  quantities  are  as  defined  in  Section  IV. 8  with  the  addition  that 
0"  =  0.2  T*  is  adequate. 

IV. 10  Modification  of  Basic  Model  by  Radial  Effects 

In  previous  sections  of  this  chapter,  it  has  been  demonstrated 
that  radial  effects  are  unlikely  to  be  either  the  principal  source  of  the 
maxpeak  effect  of  the  cause  of  smoothing  of  the  basic  maxpeak  model.  The 
main  support  is  derived  from  the  fact  that  the  maxpeak  curve  shown  by  Smith 
for  -a  core  of  radial  ratio  1.2  is  not  notably  different  from  present  curves 
for  a.  radial  ratio  of  1.68;  hence,  the  radial  ratio  is  probably  not  a 
principal  factor. 

This  viewpoint  is  corroborated  by  the  following  argument.   Let 

us  assume  that  the  basic  maxpeak  model  holds  for  an  annulus  of  core  material 

F 
and  that  H  =  -p   is  the  field  intensity  due  to  the  mmf  F  applied  to  core  path 

length  J?,  i.e.,  to  the  annulus  at  radius  JL/21T  . 

Consider  the  condition  when  FTm  =  F^T.  Then  the  time  at  which 
K  T         UP    DN 

P  domains  collapse  is  t'  =      for  each  annulus  individually.   Since 

&    d  K , 
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1.606  /FUP  _  v   ,  „    1.606  /  DN  _  x  _      ,   .   ..    . 

K  =  — - —  [—*—   -  H  )  and  Kn  =  — r —  i-n —  -  H  )  for  annulus  i.  then  for 
u     S   v  /  .    o                d     S    X  ■  o7  ' 

w   *•  1  w     1 

T 

any  i,  t'  =  —  -   a,  constant  for  all  i.   (H  and  S  are  always  assumed 
J      '      d    2  o      w 

independent  of  radius.)  Therefore,  when  F   =  F  ,    the  collapse  of  |3 
domains  should  occur  simultaneously  in  all  annulae  which  does  not,  in 
fact,  occur  as  reference  to  Fig.  k.2   will  demonstrate.  Thus,  radial 
smoothing  alone  cannot  remove  the  deterministic  fall  of  potential  which 
occurs  in  the  basic  model. 

The  point  is  that  for  one  case  (F  =  F,)  in  each  family  of 
maxpeak  effect  curves,  the  time  of  collapse  of  f3  domains  must  be  the 
same  in  all  annulae.   In  that  case,  a,  sharp  fall  of  potential  is  pre- 
dicted if  the  "basic"  majcpeak  effect  model  is  used  for  ea.ch  annulus. 
Thus,  in  one  (usual)  case,  radial  smoothing  is  completely  ineffective. 

By  extension  of  the  foregoing  argument,  it  is  seen  that  the 
predicted  maxpeak  effect  curves  for  the  total  core  must,  under  the  de- 
scribed conditions,  be  the  same  as  that  of  any  single  annulus.  There- 
fore, the  maxpeak  model  assumed  to  hold  for  an  ideal  thin  core,  whatever 
it  may  be,  will  apply  to  a  core  independent  of  its  radial  dimensions. 
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V.   THE  MAXPEAK  EFFECT  FOR  WIDER  RANGE  OF  F 

V.l  Experimental  Results 

In  the  previous  sections,  we  have  examined  the  maxpeak  effect 
for  particular  large  values  of  F .  The  only  relation  which  has  been 
found  to  hold  independent  of  F  ,  T   and  F   is  the  percentage  of  flux 
corresponding  to  maxpeak.,  namely, 


r°° 


fx     dt  =  0.83  J      xf  dt  (5.1) 


For  a  wide  range  of  F^.T,    it  is  found  that  v   -  v„  ,  x  /x_  ,  and 

DN  xp    fsp'   xp'  fsp- 

t  /t„   vary  with  F^^  as  shown  in  Fig.  5.1*   For  these  experiments, 
xp'  fsp  DN 

F  _  and  T   are  held  constant  in  order  to  sa,tisfy  the  conditions  of 
equation  ( 5.1) . 


(a)  v   -  v„  is  an  increasing  function  of  F^.T. 
v  '        xp    fsp  D  DN 

(b)  x  /x„  is  a  decreasing  function  of  F_.T. 

xp'  fsp  DN 

(c)  t  /t_  is  an  increasing  function  of  F,.,.. 

xp'  fsp  DN 


Statement  (a.)  is  in  disagreement  with  the  results  of  Smith  who, 

as  previously  noted,  finds  v   -  v_   =  constant.  The  difference  between 

xp    fsp 

Smith's  findings  and  the  present  ones  may  be  due  to  any  of  a  large  number 
of  causes  such  as  inaccuracies  in  measurement  in  either  or  both  investi- 
gations, radial  effects,  or  true  differences  in  the  magnetic  materials 
employed. 
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Statement  (b)  agrees  with  Smith  and  his  unreduced  data  appears 

to  support  statement  (c). 

V.2  Maxpeak  Effect  at  Low  F 

It  would  now  be  desirable  to  construct  a  switching  model  which 
holds  over  a  wider  range  of  F   than  the  basic  model.  For  this  objective, 
it  is  reasonable  to  attempt  to  broaden,  in  some  sense,  the  basic  model. 
Let  us  point  out  that  this  goal  has  been  only  partially  successful. 

Also,  much  of  the  following  discussion  should  be  considered 
more  speculative  than  the  preceding  work.   In  particular,  the  members  of 
the  class  of  domain  shapes  and  wall  motions  which  could  lead  to  certain 
macroscopic  results  are  clearly  vast  in  number.  Hence,  notions  corres- 
ponding to  this  class  should  be  regarded,  primarily,  as  symbolic  or 
hypothetical. 

The  argument  which  will  be  examined  is,  briefly,  that  the  large 

values  of  x  /x_   for  small  F_.T  are  strongly  related  to  the  curvature  of 
xp'  f sp  DN 

the  l/f    versus  F  (or  v  vs.  F)  curve  of  Fig.  3»3«   It  was  shown  in  Chapter 
s  p 

III  that  curvature  of  the  v  vs.  F  curve  for  full  switching  corresponds  to 

increasing  values  of  x  t  with  F.  Within  the  confines  of  the  family  of 

P  P 

curves  defined  by 

xn  =1-2  exp   [-(Kt)n]  (5.2) 

of  which  Haynes '  and  Lindsey's  equations  are  members,  variations  in  x  t 
were  shown  to  depend  only  on  n.   In  the  range  of  F  which  is  now  being 
considered  (ca.  400-700  ma.),  the  corresponding  values  of  n  were  found  to 
be  about  1.8  to  2.2.  The  meaning  to  be  ascribed  to  non-integer  n  is  now 
considered. 
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The  ca.se  n  =  2  corresponds  to  cylindrical  domain  growth, 
according  to  Lindsey's  model,  which,  neglecting  wall  annihilations,  may 
be  equivalently  stated  L  cCt  where  L  is  the  total  domain  wall  length 

t  being,  as  usual,  time  and  step  function  drive  being  assumed.   Similarly, 

2 
n  =  3  corresponds  to  ellipsoidal  domain  growth  for  which  L  oC  t  .  Therefore 

n  =  1  will  be  assumed  to  correspond  to  constant  wall  length, i.e.,  L  oC  t  =1. 

For  1  <  n  <  2,  it  will  be  assumed  that  a.  general  domain  increases 

its  wall  length  in  accordance  with  n  =  2  (that  is,  LoC  t  )  for  some  portion 

of  its  life  and  in  accordance  with  n  =  1  (that  is,  L  <£  t  )  for  the  remaining 

portion  of  its  life.  The  resulting  rate  of  flux  reversal  in  the  absence  of 

n-1 
wall  collisions  is  proportional  to  t    where  1  <  n  <  2.   For  2  <  n  <  3> 

the  reversal  process  may  be  assumed  as  a  combination  of  motions  corresponding 

to  n  =  2  and  n  =  3« 

This  approach  is  in  agreement  with  the  mathematics  of  equation 

(5.2).   It  is  therefore  possible  to  use  a,  full -switching  curve  from  that 

family  and  construct  the  ma.xpea.k  effect  curves  by  analogy  with  the  case 

n  =  2  which  is  physically  clearer.   Fig.  5-2  presents  such  surves  for  n  =  1.8, 

It  can  be  seen  that  the  maxpeak  ratio,  x  /x„   =  1.^4-U  which  is  considerably 

xp'  fsp 

larger  than  for  n  =  2  (maxpeak  ratio  =  1.33) •   Experimentally,  n  =  1.8  corre- 
sponds to  F  =  U50  ma.  (see  Table  1,  Chapter  III  and  Fig.  3.2)  for  which 

x  /x-   =  I.67. 
xp'  fsp 

The  argument,  summarized,  is  that  small  values  of  F   corresponds 
to  n  <  2  and  are  associated  with  large  maxpeak  ratios.  The  assumption 
concerning  the  meaning  of  non-integer  n  lea.ds,  indeed,  to  larger  maxpeak 
ratios  as  n  decrea.ses.  The  maxpeak  ratio  determined  from  the  model  is  not 
as  large,  however,  as  experiment  indicates. 
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The  equations  corresponding  to  this  model  are: 

(Term  l) 

.    -  (K,t_)n     K,     /t,\n-l 
n-1       d  d       d  •   f  d\ 


2(n-l)  K^"1  e     ^         -     ^    K(t) 

u     N  u' 


(Term  2) 

(5.3) 


i  *.  to"  - « *  ((■  ■  r  ■  < 


_  _d 

T 

U' 


K  T 

Term  1  applies  for  0  <  t ,  <  — - — 

-  d  -  K ,  * 
d 


K  T 
Term  2  applies  for  0  <  t.  < 


-  d  -  2  K, ' 
d 
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CONCLUSION 

1.  General 

The  principal  objective  of  the  present  investigation  was  the 
creation  of  a  quantitative  engineering  model  of  a  subclass  of  partial 
switching  effects  called,  briefly,  the  "maxpeak  effect". 

Only  scattered,  uncoordinated  experimental  data,  has  been  available 
concerning  partial  switching  effects  in  square-loop  ferrites.   The  present 
effort  is  believed  to  be  the  first  attempt  to  reduce  any  partial  switching 
effect  to  a  quantitative  engineering  model. 

In  the  pursuit  of  the  primary  goal  of  the  investigation,  some 
auxiliary  results  were  obtained  which  are,  perhaps,  of  interest.   These  are 
considered  under  the  headings  "Full  Switching"  and  "Velocity  Distribution". 

2.  Full  Switching 

The  models  of  flux  reversal  by  domain-wall  motion  as  formulated 
by  Haynes  (1958)  and  Lindsey  (1959)  have  been  used  as  the  starting  point. 
These  were  generalized  into  a.  family  of  equations  which  has  the  desirable 
property  that  the  curvature  of  the  l/f  vs.  F  curve,  which  had  been  experi- 
mentally  found,  can  be  duplicated  by  adjustment  of  the  generalizing  parameter, 
n.   It  was  shown  that  the  value  of  n  which  fits  a  given  experimental  full 
switching  curve  can  be  conveniently  calculated  from  the  quantity,  x  t  ,  which 
is  found  to  be  independent  of  H,  H  ,  and  S  . 

Since  n  can  then  be  quite  simply  and  unambiguously  derived  from 
experimental  data,  a  distinction  between  the  results  of  Haynes  and  of  Lindsey 
can  be  made.  No  clear  distinction  had  previously  been  demonstrated  because 


1.   Smith's  (i960)  collection  of  data,  is  a  notable  exception. 
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all  comparisons  were  dependent  on  the  assignment  of  the  parameters  H 

and  S  .  Thus,  "by  normalization  of  the  switching  curves  given  by  the  two 
w 

formulations,  all  practical  distinction  between  them  could  be  obliterated. 
By  the  new  method  of  comparison,  it  is  shown  that  Lindsey's  model  (n  =  2) 
corresponds  more  closely  to  experiment  in  the  region  H  <  H  <  3  H  . 

An  interesting  aspect  of  the  general  family  of  switching  equations 
is  that  the  question  of  the  physical  interpretation  (if  any)  of  non-integer 
values  of  n  arises.   By  analogy  with  the  meaning  of  n  =  2  and  3  from  Lindsey's 
and  Haynes '  models  respectively,  it  is  proposed  that  the  (non-zero)  curvature 
of  the  reciprocal  switching  time  curve  may  be  due  to  changes  in  the  average 
domain  shape  as  a  function  of  driving  field. 

This  family  of  switching  equations  has  the  additional  property 
that  all  usual  definitions  of  switching  time  are  related  by  an  easily 
computed  constant. 

3.  Maxpeak  Effect  Model  for  Limited  Range  of  F-. 

The  construction  of  a  model  of  the  maxpeak  effect  has  been  con- 
sidered. The  only  strict  invariant  in  the  maxpeak  effect  was  found  to  be 
the  flux  level  at  which  maxpeak  occurs.   It  was  argued  that  the  maxpeak 
effect  is  predominately  a.  property  of  the  magnetic  material  itself  rather 
than  the  shape  of  the  specimen  since  a  thin-walled  core,  examined  by  Smith 
(i960),  has  a  response  nearly  identical  to  a  thicker -walled  core.  A 
"basic  model"  is  constructed  by  extension  of  the  details  of  Lindsey's  model 
of  full  switching.   In  particular,  it  is  assumed  that  domain  walls  which 
are  not  annihilated  during  the  UP  drive  (|3  domains)  are  forced  to  retrace 
their  route  during  the  DOWN  drive,  and  that  renucleation  occurs  at  the 
center  from  which  the  wall  originated. 
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It  was  demonstrated  that  the  basic  ma,xpeak  model  predicts  the 
maxpeak  amplitude  ratio,  x  /x   ,  very  accurately  for  values  of  F 
(viz.,  about  2H  )  for  which  Lindsey's  equation  (n  =  2)  is  considered 
to  hold.  On  the  other  hand,  the  shape  of  the  maxpeak  effect  curves 
predicted  by  the  basic  model  is  poor  due  to  a  sharp  fall  of  amplitude 
which  occurs.  This  fall  is  traced  directly  to  the  assumption  of  Menyuk 
and  Goodenough  (1955)  (and  retained  by  Haynes  and  Lindsey)  that  the 
domain  wall  velocity  for  a  given  driving  field  is  constant  (see  Section 
k,    Conclusion). 

Therefore,  a.  new  model  was  proposed  which  incorporates  a. 
Gaussian  distribution  of  domain  wall  velocity.  The  mean  velocity  is 
assumed  to  be  the  velocity  which  would  apply  in  the  Menyuk -Goodenough 
theory.  A  standard  deviation  of  0.2  of  the  mean  time  of  collapse  of 
3  domains  has  been  found  to  be  reasonable. 

The  maxpeak  model  including  velocity  distribution  predicts 
switching  curves  whose  shape  conforms  satisfactorily  to  experimental 
results.  The  accuracy  of  the  fit  between  two  curves  is  measured  by 
the  average  absolute  value  of  the  error  in  percentage  of  peak  amplitude. 
On  this  basis,  the  error  between  the  theoretical  and  experimental,  maxpeak 
curve  is  about  6.5$>° 

It  wa.s  demonstrated  that  the  introduction  of  a.  velocity 
distribution,  as  described  above,  into  the  full  switching  equations 
results  in  only  negligible  alteration  (about  1$  in  the  switching  curve). 
Therefore,  the  full  and  partial  switching  models  are  consistent  after 
the  velocity  distribution  is  introduced. 
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Two  further  points  have  been  clarified  concerning  the  signifi- 
cance of  ra.dial  effects  within  the  toroidal  sample.  First,  the  question 
of  whether  radial  variation  of  the  driving  field  can  he  the  cause  of  the 
maxpeak  effect  has  been  investigated.   In  extensive  calculations,  (Section 
IV. 6)  a  member  of  the  family  of  full  switching  curves  (either  n  =  2,  3  or 
h)   was  assumed  to  hold  for  each  annulus  of  core  material.  The  chosen 
expression  for  flux  reversal  was  then  integrated  over  path  length  in  the 
magnetic  material  (from  inner  to  outer  core  circumference).   It  was  found 

that  no  rational  values  of  H  and  S  could  be  assigned  which  were  reasonably 

o      w 

accurate  for  both  the  full  and  partial  switching  cases.   It  was  concluded 
from  this  investigation  that  radial  effects  are  extremely  unlikely  to  be 
the  sole  cause  of  the  maxpeak  effect  and  that,  in  fact,  they  are  likely 
to  be  of  only  minor  significance. 

Secondly,  the  possible  equivalence  of  the  radial  variation  of 
the  drive  field  with  the  wall  velocity  distribution  chosen  for  the  maxpeak 
effect  was  examined.   Such  equivalence  appeared  to  be  unlikely  since,  as 
previously  noted,  Smith's  maxpeak  curves  for  a  thin  core  (OD/lD  =  1.2) 
are  not  visibly  different  from  those  of  a  thicker  core  (OD/lD  =  1.68)  used 
in  the  present  experiments .   It  is  more  probable  that  the  maxpeak  effect 
arises  from  the  material  properties,  one  of  which  is  a.  distribution  of 
domain-wall  velocities  for  a,  given  driving  field. 

h.      Maxpeak  Effect  Model  for  Wider  Range  of  F 

The  maxpeak  effect  model  discussed  to  this  point  is  applicable 
only  in  a  limited  range  of  F   .   Some  success  at  extension  of  the  basic 
maxpeak  model  to  a  wider  range  of  Fn  was  achieved  by  considering  the 
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implications  of  the  curvature  of  the  l/Z     vs .  F  curve .   This  curvature 

s 

corresponds  to  variations  in  n,  which  is,  perhaps,  related  to  the  average 

domain  shape  (see  Section  II,  Conclusion).   Since  the  wall  length,  L,  of 

an  active  domain  changes  according  to  t    for  the  full  switching  models, 

it  is  assumed  that  L  <£  t    for  non-integer  n.  Thus,  the  model  is  at 

least  mathematically  consistent  for  full  switching  and  for  the  maxpea,k 

effect  over  wide  ranges  of  F___. 

DN 

An  example  was  given  for  the  case  n  =  1.8.   It  was  found  that 
the  maxpeak  ratio  has  the  proper  direction  of  change  but  it  is  not  as 
large  as  experiment  indicates.   Nevertheless,  the  model  is,  perhaps,  of 
practical  interest  and  probably  of  some  theoretical  interest. 


5.  Velocity  Distribution 

It  had  been  assumed  in  the  Menyuk-Goodenough  theory  of  flux 
reversal  that  the  domain  wall  velocity  for  a  given  driving  field  is 
constant  (in  any  particular  specimen  of  material).   Perhaps  the  most 
interesting  auxiliary  result  of  the  present  investigation  is  that  this 
assumption  is  not  legitimate. 

The  conclusion  that  a  velocity  distribution  is  necessary  was 
reached  by  consideration  of  a  partial  switching  model.   It  is  to  be 
expected  that  a  partial  switching  model  requires  a  more  accurate 
representation  of  domain  activity  than  does  a  full  switching  model, 
the  latter  being  a  special  case  of  the  former. 

It  was  found  that  the  basic  maxpeak  model  predicts  a  large, 
deterministic  fall  of  potential  which  was  not  experimentally  observed. 
The  source  of  this  fall  was  clearly  traceable,  from  the  theoretical 
viewpoint,  to  the  formerly  cited  constant -velocity  assumption  of  Menyuk 
and  Goodenough. 
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The  possibility  that  the  fall  of  potential  was  not  experimentally 
observed  because  of  smoothing  by  radial  effects  was  considered.   It  was 
shown  that  such  smoothing  could  not  account  for  the  observations.  Addi- 
tional support  for  this  view  was  provided  by  comparison  of  a  maxpeak 
curve  which  had  been  observed  by  Smith  (i960)  for  a.  thin-walled  core 
with  present  experiments  on  a  thicker -walled  core.   Since  the  two  maxpeak 
curves  are  extremely  similar,  it  was  concluded  that  radial  effects  are 
not  responsible  for  any  substantial  smoothing. 

Therefore,  the  smoothness  of  the  maxpeak  effect  curves  is  due, 
it  was  argued,  to  a  material  property,  namely,  the  distribution  of  domain 
wall  velocity.  A  normal  distribution  having  the  mean  wall  velocity  equal 
to  the  Menyuk-Goodenough  value  and  a  standard  deviation  of  0.2  of  that 
velocity  was  assumed.  This  distribution  was  demonstrated  to  smooth  the 
maxpeak  effect  curves  in  approximately  the  desired  way.   Furthermore, 
the  introduction  of  this  same  distribution  into  the  full  switching  model 
was  shown  to  produce  only  negligible  change  in  the  predicted  full -switching 
curve. 

It  is  therefore  concluded  that  strong  evidence  has  been  provided 
for  the  existence  of  a,  distribution  of  domain  wall  velocity. 

6.  Further  Work 

The  maxpeak  models  which  have  been  proposed  are  based  on  the 
assumed  activity  of  the  so-called  p  domains.   It  would  therefore  be  of 
considerable  interest  to  experimentally  determine  whether  the  type  of 
wall  motion  hypothesized  for  these  domains  does,  indeed,  exist. 
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A  clearer  understanding  of  the  domain  activity  corresponding  to 
curvature  of  the  inverse -switching -time  curve  would  also  be  of  interest. 
This  information  would  be  invaluable  to  the  formulation  of  a  better 
model  of  switching  both  in  the  full  and  partial  switching  cases. 
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APPENDIX  I 


EXPERIMENTAL  METHOD 

1.  Drive  Pulse  Amplitude  and  Length 

Pulses  corresponding  to  F  _  and  F   were  produced  by  mercury- 
wetted-contact  relay  pulsers  in  order  to  insure  re ct angularity  of  the 
pulse  waveforms .  The  RESET  pulse,  critical  only  in  amplitude  and 
duration,  was  generated  by  a  transistor  driver. 

The  coaxial  cables  used  with  the  relay  pulsers  were  calibrated 
by  the  most  accurate  available  techniques.   Four  pieces  of  cable  having 
length  ratios  of  approximately  1:2:^:8  were  cut  from  a,  single  lot  of 
RG-62/U  cable.   Lengths  were  measured  with  a  steel  tape. 

The  velocity  of  propagation  of  each  cable  was  measured  by 
applying  a  stable  variable -frequency  oscillator  to  one  end,  the  other 
end  being  either  open  or  short  circuited.  Using  a  voltmeter  to  detect 
the  nulls  and  peaks  of  the  input  voltage  and  measuring  the  oscillator 
frequency  by  means  of  the  heterodyne  with  a  precision  100  kc/s  crystal- 
controlled  harmonic  generator,  the  velocity  of  propagation  for  the  various 
cables  was  determined. 

The  average  characteristic  impedance  of  each  was  calculated  by 
use  of  the  measured  total  capacitance  and  velocity  of  propagation. 

An  accurately  measured  non-inductive  resistor  served  as  the 

termination  of  the  pulser  output.   The  amplitude  of  the  drive  current 

produced  by  the  relay  pulser  is  then  given  by  I  = —  where 

Rt  +  Ro 

V  is  the  voltage  to  which  the  cable  is  charged,  measured 
with  a  four-decimal-digit  digital  voltmeter 


70 


R   is  the  resistance  of  the  terminating  resistor  and 

R   is  the  average  characteristic  resistance  of  the 
coaxial  cable. 


The  above  measurements  allowed  the  complete  specification  of  the  drive 
pulses  in  both  amplitude  and  duration  within  the  limitations  of  the 
following  considerations. 

Besides  the  obvious  sources  of  possible  inaccuracy  (e.g., 
accuracy  of  the  capacitance  bridge),  several  other  sources  remain. 
These  are: 

(a)  the  droop  of  the  drive  pulse  due  to  cable  losses, 

(b)  the  resistance  of  the  relay  contact, 

(c)  the  effective  shortening  of  the  drive  pulse  due  to 
the  initial  inductance  of  the  core,  the  open  wire 
passing  through  it,  and  the  relay  contact  leads, 

(d)  the  reverse  field  from  the  current  induced  into 
the  other  drive  and  sense  wires  which  pass  through 
the  core  aperture  (no  more  than  one  flux  linkage  wa,s 
used  for  any  drive  or  sense  winding), 

(e)  capacitive   coupling  of  drive  pulses  to  the  sense 
winding . 

The  effect  (e)  can  be  minimized  by  placing  the  core  at  the  ground  end  of 
the  output  drive  wire.  This  effect  is  further  minimized  by  single-turn 
windings.  The  remaining  spurious  coupling  is  of  sufficiently  short 
duration  that  observation  with  a  high  speed  oscilloscope  (Tektronix 
Type  585)  allows  unambiguous  observation  of  the  main  switching  waveform. 
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The  effects  of  (d)  can  be  estimated  by  removing  or  adding  resistance 
across  an  extra  linking  wire.  This  inaccuracy  is  usually  negligible. 
The  shortening  in  (c)  is  estimated  to  be  no  greater  than  one  millimi- 
crosecond. With  respect  to  (b),  the  d.c.  resistance  is  negligible 
while  the  variation  due  to  contact  bounce  can  be  minimized  by  adjusting 
the  relay  drive  frequency  so  that  mechanical  resonances  are  avoided. 
The  drive  pulse  droop,  (a)  is  of  some  consequence,  being  about  Vjo   per 
100  mu-s.  for  RG-62/U  cable.  The  average  amplitude  of  the  pulse  can  be 
estimated  and  used. 

The  transistor  pulser  used  for  resetting  delivers  a  pulse  of 
50  mu.s .  rise  time,  15  |is.  duration  and  amplitude  of  3a.  using  2  parallel 
type  2N1072  transistors  in  the  output.  The  amplitude  is  approximately 
2.5  times  the  maximum  drive  used  in  any  other  part  of  the  cyclic  driving 
operations . 

2.  Waveform  Measurement 

Measurements  were  made  directly  on  the  oscilloscope  display  of 
the  voltage -time  waveform  of  the  switching  core.  The  peak  amplitude  of 
the  core  output  voltage  signal  was  measured  by  shifting  the  waveform  with 
a  d.c.  voltage  equal  to  the  signal  amplitude.  The  d.c.  voltage  was  then 
measured  with  a  digital  voltmeter.   Using  Tektronix  Type  585  oscilloscope, 
which  has  a  considerably  sharper  focus  than  usual,  it  can  be  demonstrated 
that  the  maximum  error  due  to  visual  reading  of  the  beam's  vertical  posi- 
tion can  be  held  to  within  0.2  mm.  which  corresponds  to  2  mv.  for  the 
vertical  gain  normally  used  here  (100  mv./cm.). 

If,  after  shifting  the  d.c.  level  of  the  signal  for  measurement 
and  repositioning  the  beam,  the  zero  level  fell  outside  the  limits  stated 
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above  (due,  for  example,  to  drifting  of  the  oscilloscope's  vertical 
position  during  measurement),  the  measurement  was  rejected. 

Various  slowly  varying  parameters  (e.g.,  tilt  of  the  horizontal 
line,  calibration  of  the  digital  voltmeter)  were  checked  before  and  after 
each  group  of  measurements.  Temperature  was  continuously  monitored,  and, 
finally,  most  measurements  were  repeated  on  subsequent  days. 

Peaking  times  were  measured  from  the  beginning  of  the  drive 
pulse  which  is  quite  unambiguously  marked  by  the  fast  rise  time  (3-5  mu.s5) 
of  the  Type  585  oscilloscope.   Unfortunately,  however,  the  broadness  of 
the  voltage  peak  results  in  poor  accuracy  (perhaps  10$)  in  this  measurement 
even  for  relatively  fast  (and,  therefore,  narrow -peaked)  signals.   Horizontal 
calibration  of  the  scope  was  made  by  reference  to  10  or  20  mc/s.  sine  waves 
from  an  accurate,  crystal  clock. 

The  limited  accuracy  of  this  measurement  must  be  allowed  for 
when  comparisons  to  calculated  data  are  made. 

3.   Size  Measurements 

Outside  and  inside  diameters  of  cores  were  measured  by  means  of 
a  micrometer -adjusted  microscope  at  a.  magnification  corresponding  to  one 
micron  per  division  of  the  micrometer  scale.  The  average  of  two  measure- 
ments taken  along  perpendicular  lines  were  used  in  various  calculations. 

The  inner  and  outer  circumference  of  the  core  used  in  the  present 
work  were  2.238  mm.  and  3»77  mm.,  respectively. 

h.      Core  Type  Used 

A  Telemeter  Magnetic  Type  50T5  core  was  used  for  all  measure- 
ments cited  in  the  test.  This  is  a.  standard  MgMn  memory  core  apparently 
similar  to  General  Ceramics  S^  material. 
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APPENDIX  II 

RADIAL  EFFECT  AND  MISCELLANEOUS  COMPUTATIONS 

1.  General 

A  number  of  experimental  computations  were  performed  with  the 
aid  of  the  ILLIAC  during  the  course  of  the  present  investigation.  The 
various  programs  employed,  totaling  over  thirty  in  number,  can  be  cate- 
gorized as  follows: 

(1)  Compute  a  table  of  x(t)  and  x(t)  for  the  general  full 
switching  equation. 

(2)  Compute  optimum  values  of  S  and  H  (and,  optionally,  0  ) 
for  fitting  of  a.  full  switching  curve  (with  or  without 
radial  effects)  to  experimental  data.  The  minimization 
procedure  used  was  variously:  steepest  descent,  Newton's 
method,  or  a  nameless  hybrid  method. 

(3)  Compute  a  table  of  x   and  x_  as  a  function  of  t  for  given 

F  ,  T  ,  F,,  and  n  with  or  without  radial  effects. 
u7   u   d' 

(h)      Compute  x_  as  a.  function  of  t  with  x^  =  ax   +  (l-a)  x 

That  is,  x_   is  the  weighted  sum  of  two  of  the  general. 

switching  curves . 

(5)  Compute  x  with  H  or  S  as  a  specified  function  of  x. 

o     w 

2.  Radial  Effects 

Integration  of  the  switching  equations  over  radius  (or  magnetic 
path  length)  of  the  core  was  commonly  performed  in  the  following  way.  If 
some  member  of  the  family  of  switching  equations  (3»^)  is  assumed  to  des- 
cribe the  switching  of  a  differential  radial  thickness  of  core  material, 
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then  the  switching  equation  including  radial  effects  is: 


=  I 


2nt 


n-1 


'K7 


h  *a 


H 


n 


exp 


'l.606  /F 


-   H 


M> 


(A  II. 1) 


where 


and 


0  is  the  remanent  flux  (webers) 

t  is  the  time  (usee.) 

F  is  the  drive  mmf  (amperes) 
H   is  the  "threshold  field"  (amp. /meter) 
S   is  the  switching  constant  (amp,  -u.sec. /meter) 

v  is  the  emf  (volts). 


Consistent  units  are  listed  parenthetically.  The  quantity  — ~ —  (  j 


w 


HI  .,  K. 


The  constant,  I.606,  is  not  essential  to  the  formula  hut  is 
included  for  completeness.  The  constant  was  obtained,  originally,  from 
Haynes '  equation  where  it  acts  as  an  adjustment  for  the  definition  of 
switching  time. 

Equation  (A  II. l)  was  integrated  by  Simpson's  Rule  with  appro- 
priate checks  to  insure  that  the  number  of  integration  steps  was  sufficient 
to  produce  the  desired  accuracy  of  the  result. 


3.   Curve  Fitting 

Given  experimental  data  from  the  voltage -time  waveform  of  the 
oscilloscope,  it  is  frequently  desirable  to  fit  the  data  to  some  family 
of  equations.   In  the  earlier  phases  of  the  present  investigation,  the 
following  method  was  used. 
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A  number  (usually  6)  of  experimental  sets  of  values 

/v.,  t.,  F."l  were  recorded.   Starting  with  trial  values  of  H  and  S  , 
\  r  r  If  o     w' 

equation  (A  II. l)  was  solved  for  the  various  |t.,  F.>  to  obtain  v*.  The 

quantity  Q  =  Z  (v.  -  v.)  was  minimized  by  adjustment  of  H  and  S  (and, 

i 
perhaps,  0  ) .  That  is,  a  "least  square"  criterion  of  minimization  was 

used. 

The  foregoing  procedure  has  the  disadvantage  that  (a)  the  com- 
putation may  be  lengthy  and  (b)  small  errors  in  the  experimental  data  may 
strongly  affect  the  results.  To  abrogate  difficulty  (b),  it  would  be 
necessary  to  use  a  larger  number  of  experimental  points  which,  unfortunately, 
increases  the  computation  time. 

It  is  possible,  after  some  experience,  to  estimate  the  values 
of  H  and  S  which  correspond  to  some  prescribed  peak  amplitude  and  peaking 
time  at  a  given  n.  A  simple  analysis  program  is  then  used  to  verify  the 
estimated  values . 

Generally,  it  appears  from  experience  that  the  most  economical 
method  of  curve  fitting  is  by  means  of  a  large  number  of  previously  pre- 
pared graphs  of  the  various  quantities  of  interest. 

k.     Maxpeak  Effect 

An  attempt  not  previously  described,  to  simulate  the  maxpeak 
effect  was  by  linear  combinations  of  portions  of  full  switching  equations 
under  time  translation.  The  results  of  this  investigation  were  negative. 
It  is  possible  to  understand  this  by  the  following  argument. 

Consider  a  full  switching  curve  of  the  general  family  displaced 
by  time  X  as  in  Fig.  A  II. 1.  Assume  that  the  portion  of  the  curve  for 
t  >  0  is  a  maxpeak  curve.  Therefore,  the  area  under  the  curve  for  t  >  0 
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should  be  about  0.8  of  the  total  area.  Then,  the  following  equations 
must  hold. 

(l)  x    =  2n[ )    e      K,       where  DOWN  drive  is  assumed 

v  '        xp       \  n  /  d 


(2)    t    -  r  =  t 

P       xp 


(3)  Kdr=  (-m  o.8)  ' 


p     \n/ 


1/n 

v 

d 


The  above  relations  are  to  be  found  in  Section  III. 2 
Combining  (2),    and  (K); 


/n-l\  1/n 

(5)     r-  W_    - 

t 
xp 

Combining  (3),  (**■),  and  (5); 

(6)  i-  (-m  o.8)l/n  = 

Kd 

r-^)1/n 

Ka 

-  t 

xp 

Rearranging  and  using  K  from  (l); 


(7)  *  t    .  2n(-iPe'^  |Sdjy»-(^0.8)V» 

v     xp  xp       V  n  /  \  n/ 


Solution  of  (7)  for  various  n  gives  the  following  results: 
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2 

3 
k 

10 

30 


X   t 

xp  xp 


0.^03 
0.629 
0.7^6 
0.956 
1.000 


Experimentally  it  is  found  that  it   is  typically  1.375. 

xp  xp 

Therefore,  it  is  not  possible  to  simulate  a  maxpeak  curve  by  means  of 

a  single  member  of  the  family  of  (3.^0,  however  translated. 

Similarly,  it  can  be  seen  that  any  linear  combination  of 

similar  curves  would  tend  to  yield  a  value  of  x  t   which  is  too  small 

xp  xp 

compared  to  experiment.   This  is,  indeed,  the  difficulty  which  is  observed 
in  attempting  to  simulate  the  maxpeak  effect  by  linear  combinations  of 
switching  curves  having  different  n  and,  also,  by  radial  effects. 

For  both  of  these  cases,  however,  it  should  be  noted  that  each 
component  of  the  total  waveform  is  not  subject  to  the  restriction  that 
80$  of  its  total  area  lie  to  the  right  of  t  =0.   Only  the  resultant 
"maxpeak  curve"  must  be  so  restricted.  Therefore,  the  above  analysis 
does  not  strictly  apply  to  the  cases  considered. 


Figure  A  II. 1 
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